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ABSTRACT 

This book offers a concise review of quantum radar theory. Our approach is pedagogical, making 
emphasis on the physics behind the operation of a hypothetical quantum radar. We concentrate our 
discussion on the two major models proposed to date: interferometric quantum radar and quantum 
illumination. In addition, this book offers some new results, including an analytical study of quantum 
interferometry in the X-band radar region with a variety of atmospheric conditions, a derivation of 
a quantum radar equation, and a discussion of quantum radar jamming. 

This book assumes the reader is familiar with the basic principles of non-relativistic quantum 
mechanics, special relativity, and classical electrodynamics. Our discussion of quantum electrody¬ 
namics and its application to quantum radar is brief, but all the relevant equations are presented in 
the text. In addition, the reader is not required to have any specialized knowledge on classical radar 
theory. 

KEYWORDS 

radar, quantum radar, quantum information, quantum sensors, quantum electrodynam¬ 
ics, quantum optics, classical electrodynamics, radar cross sections 
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Preface 


One of the major scientific thrusts from recent years has been to try to harness quantum phenomena 
to dramatically increase the performance of a wide variety of classical information processing devices. 
In particular, it is generally accepted that quantum computers and communication systems promise 
to revolutionize our information infrastructure. 

These advances in quantum information science have had a considerable impact on the de¬ 
velopment of quantum sensors. That is, sensing devices that exploit quantum phenomena in order to 
increase their sensitivity. The interaction between quantum information science and quantum sens¬ 
ing is substantial. For example, quantum sensors can be described mathematically as noisy quantum 
channels. In addition, quantum control techniques developed in the context of quantum computation 
are useful to harness quantum sensing hardware. 

It is important to remark that, while quantum sensing is not as mature as quantum computa¬ 
tion, it offers simpler engineering challenges. Indeed, a quantum computer requires a large number 
of qubits on arbitrary superpositions and the application of a large variety of gates. On the other 
hand, quantum sensors require a small number of qubits on specific entangled states and only require 
a small set of quantum operations. As a consequence, it appears that the development of quantum 
sensors offers an enticing near term option for the practical application of the quantum information 
technologies required for the construction of a hypothetical quantum computer. 

Such is the case of quantum radar, which offers the prospect of detecting, identifying, and 
resolving RF stealth platforms and weapons systems. Preliminary results show that quantum radar 
using entangled photons can provide a quadratic increase of resolution over non-entangled photons. 
Furthermore, the “effective visibility” of certain targets is increased if observed with a quantum 
radar instead of a classical radar. Also, quantum radar is more resilient against the use of jamming 
countermeasures. And perhaps more important, a quantum sidelobe structure offers a new channel 
for the detection of RF stealth targets. 

This book offers a concise review of quantum radar theory. Our approach is pedagogical. The 
motivation for quantum radar is presented in Chapter 1. The equations that describe photons and 
their interactions with matter are discussed in Chapter 2. The reflection of photons by a mirror 
is considered in Chapter 3, which offers a novel quantum description of specular reflections using 
quantum electrodynamics and interferometry. In Chapter 4 we review classical radar and associated 
cross sections. Quantum radar theory is discussed in detail in Chapter 5, with special emphasis 
on interferometric quantum radar and quantum illumination devices. Finally, Chapter 6 presents a 
detailed discussion of quantum radar cross sections as well as examples for spherical and rectangular 
targets. This book offers some new results, including an analytical study of quantum interferometry 
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in the X-band radar region with a variety of atmospheric conditions, a derivation of a quantum radar 
equation, and a discussion of quantum radar jamming. 

This book assumes the reader is familiar with the basic principles of non-relativistic quantum 
mechanics, special relativity, and classical electrodynamics. Our discussion of quantum electrody¬ 
namics and its application to quantum radar is brief, but all the relevant equations are presented in 
the text. In addition, the reader is not required to have any specialized knowledge on classical radar 
theory. 


Marco Lanzagorta 
July 2011 




Acknowledgments 


I am grateful to Keye Martin for his advice and encouragement. I am also thankful to James Mc¬ 
Cracken, James Troupe, and Jeffrey Uhlmann for reviewing this manuscript and for providing valu¬ 
able corrections, suggestions, and insights. Special thanks to Bryan O’Gorman, who produced the 
first version of the radar cross section simulation codes. And of course, I deeply appreciate the en¬ 
couragement and perseverance provided by my editor Mike Morgan. And last but not least, I am 
thankful to Pam Lamb because her administrative support makes my life so much simpler. 

I kindly acknowledge the support received by the Office of Naval Research Quantum Infor¬ 
mation Science Basic Research Challenge program and ITT Exelis. 


Marco Lanzagorta 
October 2011 




1 


CHAPTER 1 

Introduction 


Radar is a standoff sensing technology that was originally developed during the first half of the 20th 
century [105]. While the reflection of radio waves by solid objects was first observed by Heinrich 
Hertz in 1886, it was not until the years that followed World War One that scientists realized that 
radar could become a new battlefield technology. During the inter-war period, all major military 
powers such as the US, the United Kingdom, France, Germany, Japan, Italy and the Soviet Union, 
performed radar research. 

The first pulse-based radar system, able to detect and estimate the range of a target, was 
developed at the US Naval Research Laboratory in 1934. Subsequently, radar proved to be an 
important tool during World War Two [85]. For example, the British were the first to exploit radar 
to detect airborne targets, giving them advance notice of incoming aircraft attacks. Arguably, radar 
gave a crucial advantage to the Allies during the Battle of Britain [11]. Since then, radar has become 
a widespread technology not only in the modern battlefield, but also in the civilian world. 

1.1 THE NEED FOR IMPROVED RADAR SYSTEMS 

There is a critical need to improve the performance of current radar systems and increase the 
probability of detection, identification, and discrimination of stealth platforms and other objects. In 
particular, improved standoff sensors are required for the following applications: 

Military: Stealth technology was exclusive to the US for the last decades of the 20th century [1,86]. 
However, stealth technology is becoming widespread. For example, China recently unveiled 
their plans to build their J-20 stealth fighter [100]. As a consequence, the ability to provide an 
early warning system against an enemy attack is severely diminished. 

Space Exploration: There is an incredibly large amount of natural and man-made debris surround¬ 
ing our planet [54]. From small rocks to old satellites, these objects pose a great risk to the 
manned and unmanned exploration of space. For instance, the STS-118 was damaged by a 
small object, producing a hole several centimeters in diameter in its hull [64]. 

Planetary Defense: The past few years have witnessed a renewed interest on mitigation strategies 
for planetary defense (tracking asteroids, comets, and other astronomical bodies which could 
be on a collision path with Earth) [78]. In this regard, the Shoemaker-Levy 9 event, where 
a comet impacted the surface of Jupiter, made clear the possibility of further extinction level 
events on our planet [107]. 
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As we will see, quantum information technology may be the key to improved long range 
sensor systems. 

1.2 FROM QUANTUM INFORMATION TO QUANTUM 
SENSORS 

One of the major scientific thrusts in recent years has been to try to determine the many ways in 
which quantum phenomena could be harnessed to dramatically increase the performance of a wide 
variety of classical information processing devices. In particular, it is generally accepted that quantum 
computers and communication systems promise to revolutionize our information infrastructure [81]. 

Arguably, the main motivation of most quantum information science efforts has been the 
creation of a quantum computer able to run Shor’s algorithm for crypto-analysis applications. Even 
though a full-blown quantum computer able to factorize a large number is still far into the fu¬ 
ture, many theoretical and experimental results have shown that it is possible to harness quantum 
phenomena for computation and communication purposes. 

As a side effect, advances in quantum information science have had a considerable impact on 
the development of quantum sensors', sensing devices that exploit quantum phenomena in order to 
increase their sensitivity. Examples of quantum sensors recently proposed by the scientific community 
include magnetometers, photodetectors, and gravitometers. 

Furthermore, the interaction between quantum information science and quantum sensing has 
been substantial; for example, quantum sensors can be described mathematically as noisy quantum 
channels. In addition, quantum control techniques developed in the context of quantum computation 
are useful to harness quantum sensing hardware. 

It is important to remark that while quantum sensing is not as mature as quantum compu¬ 
tation, it offers simpler engineering challenges. For instance, a quantum computer requires a large 
number of qubits on arbitrary superpositions with coherence times long enough to perform complex 
computations using a large variety of gates. On the other hand, quantum sensors require a small 
number of qubits on a specific entangled state and only require a handful of quantum operations. 

As a consequence, it appears that the development of quantum sensors offers an enticing near 
term option for the practical application of the quantum information technologies required for the 
construction of a hypothetical quantum computer. 

1.3 QUANTUM RADAR TECHNOLOGY 

Quantum radar offers the prospect of detecting, identifying, and resolving RF stealth platforms 
and weapons systems. In the electronic battlefield, quantum radar may become a revolutionary 
technology just as RF stealth technology was during the last three decades of the 20th century. The 
use of quantum radar for planetary defense and space exploration is also worth considering. 

By exploiting entanglement, quantum radar offers the prospect of enhanced target detection 
capabilities. These systems rely on quantum states of light (photons) sustained on an entangled 
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superposition. One half of these states is sent towards the target and the other remains in the 
receiver. Signal detection is enhanced by exploiting the existing correlation between the emitted 
photons bounced back by the target and the ones kept inside the radar. As such, quantum radar offer 
the possibility of detecting, identifying, and resolving stealth targets. 

Preliminary results show that quantum radar using entangled photons can provide a quadratic 
increase of resolution over non-entangled photons [67, 97]. Furthermore, the “effective visibility” 
of a target is increased if observed with a quantum radar instead of a classical radar, and a quantum 
sidelobe structure offers a new channel for the detection of RF stealth targets [61]. 

Although there are practical engineering challenges, it is fair to say that some members of 
the scientific community have shown skepticism about the prospect of building a quantum radar. 
However, to date there are no theoretical or experimental results that suggest that photons in the 9 
GHz frequency cannot be entangled, transmitted, and detected. 

Furthermore, a group at the University of Leiden in the Netherlands has already proposed a 
mechanism to generate entangled microwave photons using quantum dots [27]. This paper is cited 
in the celebrated review on quantum entanglement by the Horodecki group [47]. On the other 
hand, the Universidad del Pais Vasco in Spain has developed theoretical models for single photon 
detectors operating in the microwave regime [88, 89]. In addition, research performed at NRL has 
shown that a theoretical quantum radar operating in the 9 GHz regime is able to enhance target 
detection even if atmospheric attenuation is taken into consideration [97]. Also, a new analysis 
presented in Chapter 5 shows that supersensitive quantum interferometry is feasible up to about 300 
km even in the presence of an attenuating atmosphere. Both of these results were performed using 
the entanglement attenuation formulas derived by the Mitre group [32, 33] based on a quantization 
scheme of the electromagnetic field in an attenuating environment [52, 68]. 

It is important to note that photons in the 9 GHz frequency correspond to the radar-microwave 
regime. Indeed, the X band in radar corresponds to the 8-12 GHz region and is extensively used 
for missile guidance, marine radar, weather, ground surveillance, and airport traffic control [31, 59, 
95, 96]. Therefore, a quantum radar operating in this band would be extremely useful in a variety of 
operational environments of civilian and military interest. 

1.4 THE QUANTUM RADAR RESEARCH COMMUNITY 

There are several groups that are currently performing research on topics related to quantum radar 
and standoff quantum sensing. Among these Jeff Shapiro and Seth Lloyd at MIT [34, 48,67,93,94], 
James Smith at the US Naval Research Laboratory (NRL) [97], Jonathan Dowling at Louisiana State 
University [26, 50, 65], John Myers at Harvard University [71, 79], Yanhua Shih at the University 
of Maryland [77], Robert Boyd at the University of Rochester [72], Usha Devi at the University of 
Bristol [102], Enrique Solano at the Universidad del Pais Vasco [88, 89], Ned Allen at Lockheed 
Martin [2], James Troupe at the Naval Surface Warfare Center (NSWC) at Dahlgren, and Marco 
Lanzagorta at ITT Exelis [61, 62]. Furthermore, the Defense Advanced Research Projects Agency 
(DARPA) promoted the interest on quantum radar during their program on Quantum Sensors [17] 
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and the Office of Naval Research (ONR) recently held a technical workshop to discuss the feasibility 
of quantum radar [82]. 

1.5 SUMMARY 

Quantum radar is a promising technology that could have a strong impact on the civilian and 
military realms. Although quantum sensing technology is not as mature as quantum computation, 
the engineering of a full-blown multi-qubit quantum computer is a much harder challenge than the 
design of deployable quantum sensors. As a consequence, quantum radar is a high risk, high payoff 
proposal. By any means, the concept of a quantum radar is worth investigating in further detail. 
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CHAPTER 2 

The Photon 


The photon is the elementary excitation of the quantum electromagnetic field. That is, electromag¬ 
netic fields are made of photons. Quantum Electrodynamics (QED) is the relativistic theory that 
describes the dynamics of photons and its interactions with matter [6, 40, 90, 106]. High precision 
experiments at the European Organization for Nuclear Research (CERN) have shown that, to date, 
QED is the most accurate physical theory known to man [24]. 

In this chapter and the next we will briefly review the basic structure of quantum electromag¬ 
netic fields as they relate to non-relativistic matter fields. This specific class of applications of QED 
is usually referred as quantum optics [13,68,73,92]. A detailed discussion of QED or quantum optics 
is outside the scope of the present book, and the reader is encouraged to pursue the rich literature 
available on the topic. Our emphasis is to present those equations that are required to understand 
quantum radar mechanisms such as pulse generation and propagation, as well as the interaction 
between a quantum radar pulse and a non-relativistic material target. 


2.1 MAXWELL EQUATIONS 

The electromagnetic fields in vacuum obey Maxwell equations: 


which can be rewritten as: 


V • E 

V B 

V x B 

V x E 


0 

0 

1 9E 


dt 

1 9B 


c dt 


V 2 E- 

v 2 b- 


1 9 2 E 

c 2 dt 2 

1 9 2 B 

c 2 dt 2 


0 

0 


( 2 . 1 ) 


( 2 . 2 ) 


to explicitly describe electromagnetic waves [8, 51]. 

The electromagnetic fields can be expressed in terms of vector and scalar potentials, A and O, 
as follows: 


B 

E 


V x A 
— V<E> - 


1 9A 

c dt 


(2.3) 
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Furthermore, the potentials are arbitrary to the extent of a scalar function A: 


A —» A' = A + VA 

<F —» 


A' = A + VA 



(2.4) 


c dt 


This is known as the gaugefreedom or gauge invariance of the electromagnetic fields [51, 90]. 


2.2 ELECTROMAGNETIC QUANTUM FIELDS 


The canonical quantization of the electromagnetic field replaces the classical fields by quantum 
operators: 


(2.5) 


E 

B 


E 

B 


that satisfy Maxwell equations [68, 73, 92]. As a consequence, those relations between the electro¬ 
magnetic fields that are true in classical electrodynamics, continue to hold in quantum electrody¬ 
namics. 

For example, the reflection and refraction of a beam of light traveling across the boundary 
of two media with different electromagnetic properties can be computed using Maxwell equations 
and setting the appropriate boundary conditions for the fields [51]. Consequently, excitations of 
the quantum electromagnetic fields (i.e., photons) obey the exact same properties when they cross a 
boundary. 

In a cavity of volume V , the electromagnetic fields are made of elementary quantum excitations 
of momentum k, frequency a>k, and polarization X. The quantum fields can be expressed as Fourier 
expansions in terms of these elementary excitations as follows: 


E = T k e -^ +!k - r + h.c. 


( 2 . 6 ) 


B 



where h.c. stands for the hermitian conjugate, is the polarization vector (with X = 0, 1), and: 



(2.7) 


The annihilation and creation operators, a and 4, satisfy the commutation relations: 


[4,1,4',i'l = 0 
= 0 

[4,1 i «k',l'] = ^kk' 4l' 


( 2 . 8 ) 
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In terms of the annihilation and creation operators, the operator that corresponds to the 
electromagnetic radiation energy is given by the Hamiltonian: 

-Hr = \Jdv(e oE-E+i-BB^ (2.9) 

= I] ha >k + \ 

k,A V 

The electromagnetic quantum fields are best defined on an energy eigenvalue occupation 
number Fock space. That is, the eigenvectors are of the form: 

l«kA> (2-10) 

which represents a state of n elementary quantum field excitations with momentum k and polariza¬ 
tion X. 

The creation and annihilation operators are used to increase or decrease the number of exci¬ 
tations on a quantum state: 

l«k.i) = V«k,A. + 1 l«k,A. + l) (2-11) 

«k,A l«k,i) = l«kA - ! > 


Therefore, the creation operators can be used to define any quantum state of light from the quan¬ 
tum vacuum state |0). For example, the state of n elementary excitations with momentum k and 
polarization X is given by: 


l«k,k) 



( 2 . 12 ) 


For most applications of interest to quantum radar, the electromagnetic radiation moves as a 
beam in a straight line. For these problems, it is useful to use the continuous-mode quantization of 
the fields instead of discrete cavity modes [68, 92]. Here, the quantization volume is made infinite, 
the sum is converted into an integral: 


and the operators change to: 



2V 

(27T) 3 




(2.13) 


(2.14) 
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The one-dimensional continuous mode variable is usually taken to be co along the direction 
of propagation of the light k, and the creation and annihilation operators can be simply denoted by 
<T(&>) and a{oo). Thus, the non-zero commutation relation can be written as: 


[a(co), a^(co)] = 8(a> — co') 


(2.15) 


Furthermore, in continuous mode quantization, the equal time commutation relations for the 
actual electromagnetic fields are: 


(2.16) 


[Ei(r,t), Bj(r',t)] = -ihc 2 -&\t - r') 

01 

Bi(r',t)} = 0 

Ej(r',t)] = 0 
[Bi(r,t), = 0 

where ( i, j, l) form a cyclic permutation over the three vector components of the electromagnetic 
fields. As a consequence, orthogonal components of E and B can not be measured simultane¬ 
ously [92]. 

2.3 THE PHOTON 

In the radiation gauge, the vector and scalar potentials satisfy: 


V - A = 0 

< 1 > = 0 


(2.17) 


and the electromagnetic fields are given by: 


E = - — 
dt 

B = V x A 


(2.18) 


In this gauge, Maxwell’s equations are reduced to a single wave equation for the vector potential A: 


V 2 A - 


a 2 A 


= o 


The canonical quantization of the vector potential leads to: 

A = ^ j 


A k „ = 


(2.19) 


( 2 . 20 ) 


Vico 


where: 


( 2 . 21 ) 





with: 
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co = |k| ( 2 . 22 ) 

and the field components satisfy the orthogonality condition: 

J • A^ d^x = (2.23) 

In the above equations, e“ are the two field polarization vectors (with a = 0, 1). Because of 
the radiation gauge condition we have that: 

V • A = V • £ (A ka d ka + AlX a ) = 0 (2-24) 

k,^ 

and as a consequence: 

k • e“ = 0 Vet, k (2.25) 

This is the transversality condition which establishes that the electromagnetic fields are orthogonal 
to the direction of propagation of the electromagnetic wave. 

Also, d^a and dfy, are the electromagnetic field potential annihilation and creation operators, 
respectively, and satisfy: 

(Nka — l|«kff|Nka) = ^N ka (2.26) 

(■^ko' IMca 1 — 1) = \/~^ka 

The quantum states of the electromagnetic potential are labeled by N^ a , that is, N elementary 
quantum excitations of the electromagnetic vector potential with momentum k and polarization a. 

In the context of quantum electrodynamics, it is customary to define the photon as the ele¬ 
mentary quantum excitations of the electromagnetic vector potential [90, 106]. However, it can also 
be defined as the elementary excitations of the actual electromagnetic fields. Indeed, the functional 
representations of these quantum fields are related by the definition of electromagnetic potentials 
and by Maxwell equations in a straightforward manner [73]. 

2.4 PHOTON INTERACTIONS 

It is a well known fact that electromagnetic fields interact with charged particles. In minimal coupling 
quantum electrodynamics, photons are neutral gaugefields coupled to charged particles [90, 106] . The 
equation that describes the interaction between charged particles and photons is given by replacing 
the standard derivative terms in the non-interacting dynamical equation of the free particle with the 
covariant derivative [90,106]. 

For instance, the Dirac Lagrangian for a non-interacting fermion of mass m described by the 
Dirac spinor fy is: 


C = — mifnfr 


(2.27) 
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where y^ are the Dirac matrices and 3^ is the standard 4-derivative: 


8 “ = l5? v 


(2.28) 


Then, the electromagnetic interaction is introduced by replacing 3^ with the covariant derivative 


D T 


— 3 (i i e A 11 


(2.29) 


where e the electric charge of the fermion and A^ is the Lorentz covariant electromagnetic potential: 

A„ = (</>, -A) (2.30) 

Therefore, the interacting fermion is described by the Lagrangian: 

C = ixjsy^D^xlr — mxlrxjs 


= ixfry^d^xlr — mxlsxjs + ex/sy 11 

That is, the interaction between the charged fermion and a photon is given by the term: 

If, Y = ex/sy 11 A^ifr 

which involves two fermionic and one photonic fields. 


(2.31) 


(2.32) 




Photon Absorption 


Photon Emission 


Figure 2.1: Minimal coupling quantum electrodynamics: the interaction between light and matter is 
described solely in terms of photon absorption or emission processes. 

The physical interpretation of this interaction term is shown in Figure 2.1: the interaction 
between a photon and a fermion is simply described as a process of photonic absorption or as a 
process of photonic emission [16, 40]. Notice that each of these Feysiman Diagrams involves two 
fermionic lines and one photonic line, in accordance with the interaction term I f,y. Therefore, all 
electromagnetic interactions boil down to processes of photonic absorption and/or emission 1 . 

^ Even though we derived relativistic expressions, these results are also valid in the non-relativistic regime. 
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The photon is well described by the equations that emerge from the quantum field treatment 
of the electromagnetic fields. However, sometimes it is difficult to conceptually understand the 
dynamics of physical phenomena involving photons. One of the reasons for this shortcoming is that, 
in strong contrast to other relativistic particles such as the electron, the photon cannot be localized. 
That is, it is mathematically impossible to build a continuity equation using localization probability 
distributions that satisfy Einstein’s special relativity (see Section 2.6). 

Furthermore, the theoretical scope of QED does not include the explicit space-time de¬ 
pendence of the physical processes that it describes [6] . For instance, QED does not describe the 
continuous evolution in time of a photon being absorbed or emitted by a charged particle. Instead, 
the theory provides predictions of dynamic observables (energy, momentum, polarization) based on 
the initial state of the particles before the interaction takes place. Therefore, QED merely provides 
the transition amplitudes between specified initial and final states of a system of particles [6, 90,106]. 

These observations enhance the challenge of properly understanding how an incident photon 
is reflected by a mirror (clearly, this is a problem closely related to the interaction between a quantum 
radar pulse and a target). Indeed, the reflection process can be reduced to a scattering problem: the 
incident photon is being scattered by one or more of the atoms that make the mirror. That is, the 
incident photon is absorbed by one or more atoms, and subsequently emitted. However, how many 
atoms scatter the photon? If we cannot localize the photon, then, what atoms are involved in the 
scattering process? Also, if specific paths make no sense in the context of QED, then, how can we 
expect to obtain the law of reflections (incident angle equal to outgoing angle of the light beam)? 
Furthermore, if the scattering process in QED does not involve any geometric angles, how does the 
scattering atom “know” that it has to send the outgoing photon on a precise direction? 

Clearly, picturing the photon as a ball scattered by the atoms in a game of billiards will not lead 
to the correct answers to these questions. As such, understanding the QED description of specular 
reflections shows how the commonly held interpretation of the photon as an electron-like particle is 
completely erroneous and misleading. In the next chapter we will discuss a formal QED description 
of specular reflections. 

2.5 THE CLASSICAL LIMIT 

Quantum electrodynamics offers the only theoretical description of electromagnetic fields that is 
fully consistent with quantum theory. It is also known that the physical properties of a quantum 
system are equivalent to those of its corresponding classical system when the relevant quantum 
numbers are large " [6, 60]. In the case at hand, classical electrodynamics is recovered in the limit of 
a large number of photons. 


^This statement is not exactly true: one could envision a maximally entangled state made of a large number of photons N. Clearly, 
the limit when N is large does not lead to a description of decoherence. Indeed, the state remains entangled regardless of the 
value of N. As a consequence, the formal transition to the classical limit needs to involve both a large number of photons and a 
decoherence process. 
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However, an infinitely large number of photons also implies an infinitely large energy Clearly, 
a more rigorous condition is required to determine when a system has to be described with quantum 
or classical electrodynamics. 

To this end, let us consider a time-dependent quasi-classical electric field E averaged over a 
short time interval At. We assume that the field has a Fourier representation of the form: 

E(x, 0 = f E(k, a))e~ i{mt - k l) dk 3 dao (2.33) 


If E is the average field over At, then the only Fourier components that will contribute to E are 
those with frequencies: 

o> < — (2.34) 

At 

as for the other components the oscillatory factor e~ lait is very small when averaged over At. In 
addition, in a unitary volume, the number of oscillators with frequency between zero and co is: 



1 

(cAf ) 3 


(2.35) 


We also know that the field energy density is proportional to E 2 . Therefore, if N m is the average 
number of photons we have that: 


N w x hu> x ( —) ~ E 2 

(2.36) 

and as a consequence: 


- E 2 c 3 

(2.37) 

Then, if the number of photons is large: 


- 7 hco 4 

1 =4- E — r- 

r J 

(2.38) 


so is the value of E 2 , but the electric field is large only in reference to the averaging time interval At. 
That is: 

|E| » (2.39) 

(cAr) 2 

This is the meaningful condition that establishes when an electric field can be treated classically 
(and a similar condition can be deduced for the magnetic field) [6]. 

For example, for experiments where At is very small, the fields are required to be very strong 
if they are to be considered classically. In practice, however, At is related to the time during which 
the field changes in an appreciable manner. Thus, for a static field At —> oo and the classicality 
condition can be satisfied by arbitrary fields. In other words, all static fields are classical fields. 

On the other hand, let us assume now a time-dependent field that changes in an appreciable 
manner over a time interval At. Then, if the field is sufficiently weak, it can never be considered as 
a classical field and has to be described using quantum electrodynamics. 
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Non-relativistic quantum mechanics prescribes that the localization of a subatomic particle (i.e., 
an electron) described by the wave function T(x) is probabilistically distributed in space [40, 60]. 
That is, we can define a probability density for localizing the particle in the volume element d}x as 
follows: 


p(x) = |T(x)| 2 

Together with the probability flux density: 

j(x) = - A- (T*(x)VT(x) - T(x)V4>*(x)) 
they satisfy the continuity equation for particle localization: 

3/0 (x) 


at 


+ V • j(x) = 0 


(2.40) 


(2.41) 


(2.42) 


A similar relation emerges for relativistic spin-1/2 particles described by the Dirac equation [40]: 

j M (x) ee -e'P(x)y /i 'P(x) = (p,j) =>■ dfij 11 (x) = 0 (2.43) 

In the case of photons, however, it is mathematically impossible to build a current J^ that 
satisfies a gauge invariant and Lorentz covariant continuity equation of the form [6]: 

3^00 = 0 (2.44) 

with a positive definite time-component (i.e., a positive definite localization probability density). 
Indeed, gauge invariance requires that the 4-vector electromagnetic field A ^ can only appear as the 
antisymmetric tensor: 

F^v = 3 ,iA v - d v A m (2.45) 

= i(k/jA v kyAfP) 

Therefore, the photon localization current has to be a positive definite function of the form: 


(2.46) 

However, because of the transversality condition 1 : 

k a F a p = 0 (2.47) 

it is easy to verify that all possible combinations of F^ v , F* v , and k^ will either cancel to zero or 
depend explicitly on k^ (which violates the positive definite condition). 

Furthermore, the fact that the photon cannot be localized also implies that we cannot define 
a consistent wave function for the photon. Indeed, the wave function is defined to characterize the 
spatial localization of quantum particles such as electrons. 

Ah is condition establishes that in an electromagnetic wave, the electromagnetic fields E and B are perpendicular to the propagation 
vector k. 
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2.7 PHOTON DETECTION 

The quantum electric field can be separated into its positive and negative frequency terms as: 

E(r, t) = E (+) (r , t) + E {_) (r , t ) (2.48) 


where: 


E (+) (r, t) 


E (_) (M) 


^ ' f'kf.( 1 ■ 

k,c* 

i J2 e“ e- i(at - k r) a ka 

k,a 

E E k 

k,oi 

-i J2 e“* e d^-k-r) 

k.o' 


(2.49) 


Thus, E (+, (r, t ) and E^ '(r, t ) involve the annihilation and creation of photons, respectively. 

A detector at the position (r, t ) measures a photon by destroying it (using, for instance, the 
photoelectric effect). Therefore, E (+) (r, t) is the only term of the quantum electric field operator 
that contributes to this measurement process. Furthermore, the transition amplitude of the detector 
signaling a photon is given by: 

o/(r, t) = </|E (+) (r, t)\i) (2.50) 

where the electric field went from the initial state | i) to the final state \f) through the annihilation 
of a photon. The transition probability is then given by: 


w f (t, t) = |fl/(r, Ol 2 


</|E (+) (r,t)b> 


(2-51) 


The final state of the electromagnetic field is of no interest to the measurement process, and therefore 
we need to sum over all the possible final states: 


W(r,t) = J>/(r,t) (2-52) 

/ 

= ^|</| E (+) ( r , t )|/)| 2 
/ 

= ^<i|E ( ->(r,0l/M/|E (+) (r,0b> 

/ 

= (i|E <-) (r, t) ■ E (+) (r, t)\i) 

= (E (_) (r, t) ■ E (+) (r, t)) 
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where we have used the fact that the final states form a complete basis of the Hilbert space: 

^l/x/l = l (2.53) 

/ 

Therefore, we can define the expectation value of the measured intensity as: 

(7(r,t)> = W{r,t) (2.54) 

= (yilE^Cr.O-E^Cr.Olyi) 

= Tr ^Pi E ( ~ } (r, t) ■ E (+) (r, f)) 

where the expectation value of the electric quantum field operator is taken on the initial state of the 
photon field |y, ) [92], and Tr is the matrix trace operator. 

Finally, it is necessary to mention that the First Order Correlation Fmiction is defined as: 

G (1) (ri,r 2 ; U - h) = (E (_) (ri, t\) ■ E <+) (r 2 , t 2 )) (2.55) 

Then, it terms of this function the expectation value of the measured intensity is given by: 

(7(r, t)) = G (1) (r, r; 0) (2.56) 


2.8 THE PHOTON WAVE FUNCTION 


We can presume that the photon in question was produced by some form of atomic decay. In such 
a case we can use the Weisskopf-Wigner theory of spontaneous emission of photons between two 
atomic energy levels [13, 92]. Working in the rotating wave approximation, one can approximate 
the state of the emitted photon as: 


iKo) = 


a —i kro 


Gk - Cl>) + i T/2 


Ik) 


(2.57) 


where co is the frequency difference between the atomic energy levels E u and 71),, ro is the position 
of the atom, T is the inverse of the lifetime of the excited state of the atom: 


and: 


with: 


F = 


Igki 2 


I = 1 4ftj3 lMl lb 

x 4?reo 3 he 3 

V]^ $ ^ 

= 2 he 0 V ' Ckl 

flab = (a\Jl\b) 


(2.58) 

(2.59) 

(2.60) 


the transition amplitude of the electric dipole moment of the atom between the states a and b. 
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Within the theoretical framework of the Weisskopf-Wigner model is possible to give an 
“operational” definition of a photon wavefunction [73, 92]. Indeed, suppose we have a single photon 
produced by atomic spontaneous emission. For added simplicity, also assume that the electric field 
is polarized in the direction ?. Then, the first order correlation function is given by: 


G (n (r,r;0) 


(ml E ( } (r, t) • E ( + } (r, t) |y 0 ) 

(yo|£ (_) (r, t)e * ■ E i+) (r, t)e \yf 
{y 0 \E ( -\r,t)E( +) (r,t)\y 0 ) 

(y 0 \E ( ~\rJ) \Yi)(Yi \) E {+) (r, t)\y 0 ) 


(2-61) 


where we have assumed that the photon states form a complete basis of the Hilbert space. However, 
because yo represents a single photon state, only the vacuum contributes to this term: 


G (1) (r,r;0) = (y 0 |F M (r,O|0)(0|F (+) (r, t)\y Q ) 

= <0|F (+) (r,0 |ko) 2 


(2.62) 


Furthermore, as discussed in the previous section, G (1) (r, r; 0) is the probability of measuring a 
photon at the detector site. Therefore, the function: 

V y (r,t) = <0|£ (+) (r,t) |ko) (2-63) 

can be understood as a photon wave function 4 . 

If we write the quantum electric field as: 

E (+) (r, t) = iJ2 e“ e - i(cot - k T) a ka (2.64) 

k,oi 


and transform the summation into an integral, we get the spherical wave function: 


where: 


Tn 


¥ y (r,0=^e|f- — 


A r 


g -(/ffl+r/2)(f-Ar/c) 


So 

A r 


CO 2 \fl\ab sin 

4tt€o c 2 A r 
|r - r 0 | 


(2.65) 


( 2 . 66 ) 


and ?7 is the angle between the electric dipole jl and Ar = r — ro [92]. Notice that To is independent 
of r in the far field. Indeed, \jl\ a b depends linearly on r, canceling the dependence on r introduced 
in the denominator of the expression for To. 

4 Of course, this statement is approximately true only within the context of the specific experiments discussed in this text in- 
volving quasi-monochromatic fields. Any attempt to generalize the concept of a photon wave function will ultimately lead to 
contradictions [73]. 
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2.9 PROPAGATION IN ATTENUATING MEDIA 

Any realistic description of a classical or quantum radar needs to account for the attenuation losses 
due to the atmosphere. We consider two major contributions to atmospheric attenuation: absorption 
and scattering [43]. Absorption is the process when the light is absorbed by the medium and its 
energy is dissipated as heat. On the other hand, scattering is the process when light is absorbed by 
the medium and re-emitted in some direction with minimal thermodynamic energy dissipation. In 
this section we will briefly discuss the description of attenuation in the classical and quantum realms. 

2.9.1 ATTENUATION OF CLASSICAL LIGHT 

In the context of classical electrodynamics, the energy flux of a parallel beam of light <t>(r) is 
attenuated due to absorption as: 

cp(r) = <J>(0) e~ Xar (2.67) 

where Xa is the absorption attenuation coefficient. If scattering losses are present, then the total 
attenuation coefficient becomes: 

Xc = Xa + Xb (2.68) 

where Xb is the scattering attenuation coefficient [91]. 

In the classical context, the process of absorption is well described through Maxwell equations 
in material media: 


V D 

= 47 rp 


V B 

= 0 


V x H 

47r 

1 3D 

= — J + 

— 


C 1 3B 

c 3 1 

V x E 

c 3 1 



along with the constitutive equations: 

D = E + 4ttP 
H = B - 4ttM 

If we assume a linear response medium, then we can write: 

P = XeE 


and therefore: 


D = eE 


€ = 1 +4nxe 


(2.69) 

(2.70) 

(2.71) 

(2.72) 


(2.73) 

(2.74) 


(2.75) 

(2.76) 

(2.77) 


with: 
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And similarly: 

B = /rH (2.78) 

where Xe is the electric susceptibility, € is the dielectric constant, and /t the magnetic permeability 
of the medium [51]. 

Using a simple molecular model based on damped harmonic oscillators, the dynamics of a 
molecular electron in the presence of an electric field E are described by: 

m ^3c + yx + co^x^j = —eE (2.79) 

where m, e, and coo are the electron’s mass, charge, and binding energy, respectively, and y is the 
damping coefficient [51]. 

Using this model for an arbitrarily large number of multi-electron molecules is possible to 
derive the dielectric constant as: 


e(<w) — 1 + 


4 jtNe 2 


E 


fj 


a)j — co 1 — icoy 


(2.80) 


where N is the number of molecules per volume, the summation is done over the number of electrons 
per molecule, coj is the frequency of the harmonic oscillator tied to the j th -electron, yj is the damping 
constant of the j th -electron, co is the frequency of the light traversing the medium, and /; is the 
degeneracy of the electronic states and satisfies: 


E^ = Z < 2 - 81 ) 


where Z is the number of electrons per molecule. 

Even though this expression for the dielectric constant is accurate, we need a quantum- 
mechanical model to describe the dynamics of the absorbing molecules. This is necessary to accurately 
determine yj, coj, and fj. Then, using Maxwell’s equations we find that the absorption coefficient 
is given by: 


2co 

Xa oc —3 
c 


vN 


AnNe 2 

> +-E 


fj 


m 


co 1 - — co 1 — icoy 
i J ' 


(2.82) 


where 3 represents the imaginary part of the argument. 

The estimation of the scattering coefficient Xb is much more difficult to do following ele¬ 
mentary principles [91]. One is required to use the Rayleigh scattering or Mie scattering formalism, 
which leads to accurate but rather complex formulae [8, 51]. Furthermore, the total scattering coef¬ 
ficient will depend on the size, density, and distribution of the scattering centers. In most cases, one 
is required to search the literature for empirical values of Xb [91]. 
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2.9.2 ATTENUATION OF QUANTUM LIGHT 

The quantization of the electromagnetic field in attenuating media is much more complicated than 
in vacuum [21, 49, 55, 69]. The consistent theoretical treatment of attenuating media requires the 
quantization of interacting electromagnetic, matter, and reservoir fields [49] . 

A full treatment of this problem is outside the scope of the present book, but we can say 
that the expressions for the quantized fields are very long and complicated. While this quantization 
procedure is so robust that it allows the prediction of the dielectric constant of the medium, the 
resulting expressions are too difficult for any analytical study, and it is very difficult to integrate 
phenomenological values of the attenuation constant into the theory. 



Figure 2.2: The attenuation of light (scattering + absorption) can be modeled using a chain of beam¬ 
splitters [52, 68]. 


However, it is possible to obtain simplified equations if one assumes that the medium is 
adequately described by a chain of beam splitters as shown in Figure 2.2 [52, 68]. The incident light 
in each beam splitter is represented by the operator d;„. The light that is scattered or absorbed by 
the medium is represented by a s . Also, the possible contributions from the medium to the light field 
are represented by b. Finally, the net output of light is described by a ou t ■ Then, the traversal of the 
chain of beam splitters modifies the annihilation quantum operator into: 

a{co) e ikz-XcUo)z/2 a(aj) + ; v / Xc(w) ( dx e (‘k-Xc(xo)l2)(z-x) £ (<u) ( 2. 83 ) 

Jo 


where: 


CO T]{(0) 

k = - 


c 


(2.84) 
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and r](co) is the refraction index of the medium [52, 68]. As usual, the expression for the creation 
operator is the Hermitian-conjugate of a(a>). 

This quantization scheme offers a simpler phenomenological theory where it is relatively easy 
to plug the experimental values for the attenuation coefficient Xc ■ However, because of the presence 
of attenuation due to an external field b(a>), the quantum light field has to be treated as an open 
quantum system [74]. 

As an example, let us consider a NOON state: 


| NOON) 


— (|A)i|0)2+|0)i|iV)2) 

V2 

\a\) N 0 A) N N 

1 -|0>i|0>2 + -^|0)i|0>2 



s/m 

(4) jv |0)i|0) 2 + 

(aj)^ (al) N 

s/m s/m. 


Vm. 

l 



(4)^10)110)2 


(2.85) 


After propagating each component across distances L \ and Li over two different media characterized 
by refraction indices rj i and p 2 > and attenuation coefficients Xcl and Xc2, we get: 


| NOON)' 


c -ia>m/c-Xr\(u>')/2)(NLi) 

s/2m. 

i 1 - ia>r)2/c-xMco)/2)(,NL2 ) 

s/2m. 

+ 1^) 


(aJ) A? |0) 1 |0) 2 

( 4 )^ 10 ) 110)2 


( 2 . 86 ) 


where |<h) is a state that represents those states that have been scattered outside the NOON basis 
due to the b(co) contribution to a{u>). 


2.10 SUMMARY 

The photon is the elementary excitation of the quantum electromagnetic field. Furthermore, the 
photon is a non-localizable particle; that is, it is mathematically impossible to build a continuity 
equation for the photon current. As a consequence, it is impossible to write down a wave equation 
for the photon. 

In spite of these complexities, it is possible to write down expressions for the probability of 
photon detection as well as an approximate photon wave function for monochromatic light. 

Finally, it is possible to build a quantum theory for the attenuation of light considering a chain 
of beam splitters. This quantization scheme produces equations where one can insert phenomeno¬ 
logical values of the attenuation coefficient to analyze the state of attenuated photons. 
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CHAPTER 3 

Photon Scattering 


As we discussed in the previous chapter, in the large photon number limit, Quantum Electrodynamics 
(QED) leads to Classical Electrodynamics (CED). Indeed, the electromagnetic quantum fields E 
and B satisfy the same Maxwell equations as the classical fields E and B [90,106]. As a consequence, 
those physical phenomena that are described by CED can also be described with QED. 

Such is the case of specular reflections, that is, the reflection of light by a mirror. By imposing 
the right boundary conditions on the classical electromagnetic fields one can use Maxwell equations 
to derive the well known “Law of Reflections” (the incident and outgoing rays of light form equal 
angles with respect to the normal at the mirror surface) [8, 51]. Therefore, it is to be expected that 
photons, the elementary excitations of the quantum electromagnetic field, will be reflected by a 
mirror following a similar law. 

Nevertheless, it is highly valuable to develop a description of specular reflections exclusively 
using QED fundamental processes of absorption and emission of photons [62]. First of all, it is an 
enlightening educational exercise to use the tools of Quantum Field Theory (QFT) and QED in a 
different context to their traditional use in high energy physics. Second, this problem makes evident 
the conceptual problems of interpreting the photon as an electron-like quantum particle. And equally 
important, the derived formulas found in this effort can be directly applied to quantum sensing. In 
particular, as we will see in the following chapters, a QED description of specular reflections is 
of fundamental importance for the understanding of Quantum Radar and Quantum Radar Cross 
Sections (QRCS) [61]. 

Unfortunately, a quantum description of specular reflections that explicitly uses fundamental 
QED processes is not readily available in the scientific literature. Indeed, it is not found in traditional 
textbooks on QFT, QED, or quantum optics, and neither in the scientific literature from recent years. 
A few classical optics books talk about the subject, but completely avoid all technical details [43]. 
And even though Richard Feynman provides an outstanding layman-level discussion to the subject 
matter in one of his books, he avoids fundamental QED processes in favor of a quantum path integral 
approach [30]. 

In order to overcome this deficiency, this chapter will provide a quantum description of specular 
reflections that exclusively uses fundamental QED processes of absorption and emission of photons . 
Our approach is pedagogical in nature and attempts to bridge the gap between our understanding of 
classical and quantum physics. Clearly, this is a topic that deserves further discussion and analysis. 


Uhis chapter is based on a paper written by the author [ 62 ]. 
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3.1 PHYSICAL PROPERTIES OF THE SPECULAR 
REFLECTION OF LIGHT 

The reflection of light by a mirror is a well understood physical phenomena [8, 43]. The Law of 
Reflections relates the incident and reflected angles of a beam of light and was first derived by Hero 
of Alexandria in 10-70 AD. The formal mathematical description of this phenomena was first given 
by Fresnel in the early 1800s using wave optics. Subsequently it was found that Fresnel’s equations 
could be easily derived using Maxwell equations by imposing the adequate boundary conditions on 
the mirror’s surface. 

However, the Law of Reflections is not the only physical property that characterizes specular 
reflections. The reflection of light by a mirror is characterized by the following five fundamental 
physical properties [8, 43]: 



Figure 3.1: The geometry of light reflected by a mirror. 

PI: The Law of Reflection. The angle of incidence 6j is equal to the angle of reflection 0 r (see 
Figure 3.1). 

P2: Coplanarity. The incident ray i, the reflected ray r, and the normal to the mirror n are found 
in the same plane n (see Figure 3.1). 

P3: Fermat’s Principle. The light travels on the path of least traversal time. That is, light travels 
over the shortest path. 

P4: In/Out Frequency Invariance. The frequency of the incident light co, is equal to the frequency 
of the reflected light co r . 

P5: Quartic Frequency Cross Section. The scattering of light by a mirror is due to Rayleigh 
scattering processes with a cross section do oc oj 4 . 

Therefore, the ultimate goal of this chapter is to adequately describe all of the above five 
physical properties in terms of QED fundamental processes of absorption and emission of photons. 
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3.2 ATOM-PHOTON SCATTERING 

The physical process of interest for the present discussion is the atomic scattering of radiation of the 
form: 

At + Yi —» A f + Y f (3.1) 

That is, an atom A,- and a photon Yi interact with each other in a process that ultimately leads to 
the atom A f and a photon yy in presumably different quantum states. However, conservation of 
energy dictates that: 

Ej + hcoi = Ef + ha>f (3.2) 

where E a refers to atomic energies and frcob to photon energies. 

Depending on the energy of the final states, there are four types of atom-photon scattering 
processes [16]: 

• Rayleigh Scattering: Low-energy elastic scattering (cot = co/). 

• Raman Scattering: Low-energy inelastic scattering (ftj, ^u>f). 

• Thomson Scattering: High-energy elastic scattering (oj, = o>j ). 

• Compton Scattering: High-energy inelastic scattering (<o/ 7 ^ ioj ). 

In the present chapter we will assume the low energy regime. That is, we will only consider the case 
of low energy photons 2 . We still need to show, however, that Rayleigh Scattering is the process that 
characterizes specular reflections. 

3.2.1 THE TRANSITION AMPLITUDE 

At the lowest order in non-relativistic perturbation theory we have two possible time-ordered dia¬ 
grams that contribute to the process of photon scattering by an atom (Figure 3.2) 5 : 

I The atom A,- absorbs the photon y,- transforming into the intermediate A n state. Subsequently, 

the atom emits a photon yy and goes to the A y state. 

Ai + Yi —* A n —> Af + Yf (3-3) 

II The atom A,- emits the photon Yf and transforms into the intermediate state A„. Subsequently, 

the atom absorbs the photon y,- and enters into the A y state. 

A{ + Yi — * A n + Yf + Yi —* Af + yf (3-4) 


- [ n this context, the energy of a photon is considered low or high relative to the ionization energy of the interacting atom. For 
quantum radar applications, we are indeed below the ionization limit. 

3 Note, however, that in relativistic quantum field theory both time-ordered diagrams would be represented by a single Feynman 
diagram. It is also worth noticing that while the photon is always a relativistic particle, we have chosen the atom to be in non- 
relativistic kinematic states. As all the particle propagators in the relevant diagrams refer to atoms, we can use non-relativistic 
perturbation theory to describe the dynamics of the atom [40]. 
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t 


Figure 3.2: The two lowest order time-ordered diagrams in non-relativistic perturbation theory that 
contribute to the process of photon scattering by an atom. 

The vertices of the diagrams shown in Figure 3.2 represent fundamental QED processes of 
absorption or emission of photons as described in Chapter 2. These vertices are described by those 
energy terms in the Hamiltonian that couple radiation and matter. From left to right, as seen in 
Figure 3.2, the four vertices in these diagrams represent: 

Vf : Aj absorbs a photon a>j and transforms into A n . 

Vj n : A n emits a photon u>f and transforms into A /. 

yf: Aj emits a photon u>f and transforms into A„. 

Vj n - A n absorbs a photon tu,- and transforms into A /. 

Explicitly, in terms of atomic and photonic quantum states, these vertices are given by: 


v a 

m 

= (A„| ® (o y" lA/)®!®,-) 

(3.5) 

Vi 

fn 

= (A/| (g) (u>f\ V e A,;) <2> 0} 


V e - 

m 

= (A n \ (g) (co f \ V e A,-) ® 0) 


v a f 

fn 

= (A/| <g) (0 V a |A„> O |o)/> 



where V a and V e are quantum operators that represent the atomic absorption and emission of a 
photon, respectively. And also, to simplify notation we have rewritten the initial state of the photon 
as: 

M = |l k „ ff/ ) (3.6) 

with coj = |k, |. 

Therefore, non-relativistic perturbation theory prescribes the transition amplitude Vfi for the 
atom-photon scattering process as follows [6, 60]: 

v - ' / V fn V ni VfnK \ 

l r \y £i _ s i + Si _ £ u) 


(3.7) 
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where is the initial energy, £is the intermediate energy for diagram I, and SjJ is the intermediate 
energy for diagram II 1 . These energies are clearly given by: 


£i = E, + hu>i (3.8) 

f 1 — F 

£' n ' = E n + h(Oi+hcof 

In minimal coupling QED, the interaction between light and matter is given by a term in 
the energy that couples a charged particle field with the quantum electromagnetic field [90, 106]. 
However, in the present case, mirrors are made of neutral atoms. The energy terms that describe 
how neutral atomic quantum fields interact with the electromagnetic quantum field involve the 
atom’s electric and magnetic multipole moments [13, 73, 92]. For the case of small frequencies, 
all transitions can be considered in the electric dipole approximation [6, 13, 73, 92]. Then, higher 
atomic electric and magnetic multipole moments can be safely ignored in the present discussion. 

The electric dipole ji is defined as: 


4 = ^ qpi (3-9) 

where qt represents the electric charge in the system with position r,- [51]. Clearly, /2 depends on the 
physical structure of the atom and may result in a very complex term. However, in our discussion it 
is not necessary to assume any specific form for jl. 

As discussed above, the only energy term that couples the atom quantum field with the 
electromagnetic quantum field is given by the electric dipole coupling: 

V = -£ ■ E (3.10) 


where fl is the atom’s electric dipole operator and E is the electromagnetic quantum field. 

The electric dipole interaction does not distinguish between absorption and emission pro¬ 
cesses. Thus, the emission and absorption couplings are exactly the same: 


V a = V e = V 

In the radiation gauge we can write E as: 



— ^ (Eka«k(/ + Ek(/'k« ) 

k.a 


(3.11) 

(3.12) 


where: 


Eka = iM Ak„ (3.13) 

= iyfl^o e“ 

^It can be shown that the energy terms in the denominator of the transition amplitude form the time-ordered description of the 
atom propagator in relativistic QED [40]. 
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and therefore: 

V = -jl ■ E (3.14) 

= — M ^ (^ko^ka + Ek a ^kff) 

k,a 

= sj2 ™ ( e “ e~ i(at - kl) a kci - e“* e i{M, ~ kt) aQ 

k,a 

= -i (A • e“ e- i(a>t - k ' c) 4 a - A ■ e“* e' ( “^ k ' r) a*,) 

k,a 

Notice that in the above expression, the operator jl only acts on atomic states, while the operator 
4 q. only acts on photon states. As a consequence, they commute with each other. 

For added simplicity, we can consider the case where the photon states are described by plane 
waves and the wavelength is larger than the size of the scattering atoms. This is the electric dipole 
approximation and implies that: 



, t r lit 

e k 1 + i — r + ... Rs 1 

A, 

(3.15) 

where |r is consider to be of the order of the size of an atom [6, 73]. 

Therefore, the first transition amplitude of relevance is given by: 


y a 

m 

= (A„|®(0| V |A,-)®|ffli) 

= ~i ^ V27 rw ( A n \Jl\Ai ) ■ e“ e~ lmt (0|4J<»i) 

k,a 

+i Y]^ 2nco ( A n\it\Ai) ■ e“* e' M (0|d^[u>;) 
k,a 

= -iyfTjtati (A„|/t|A,) ■ e'’ e~ i(0it 
= —iy/lizcoi (e' • finij e~ la>it 

(3.16) 

where we have defined: 

/tab — (A fl | ft 1 Aft) 

and A a , Ab represent atomic states. Similarly, for Vj n we get: 

(3.17) 


V fn = ( e/ * ' F«/J e '“ ft 

(3.18) 


And similar expressions can be derived for VT and V° n . 


3.2.2 THE SCATTERING CROSS SECTION 

In terms of transition amplitudes, the scattering cross section is given by [60]: 

_ or f d£l 
da = 2n\V f i\ 2 ——y 


(3.19) 
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and therefore, for the atom-photon scattering case in consideration we have: 

e f *)ti*ni ■ e‘) , (/!/„ ■ e')(/!,„■ • e-^*) 


flf(7 = 


| (A/>; 

n l 


to)ni to)i 


■e‘) ■ e')(/t„ 

-1--- 

U) n f T COi 


2 3 

COjCOj 

h 2 c 4 


dC2 


(3.20) 


where we have defined: 


_ E a - E b 

to)ab — » 


(3.21) 


The above equation reveals a few interesting facts. First of all, the cross section is zero if 
the electric dipole is orthogonal to the polarization vector of the photon. In other words, the atom 
cannot absorb or emit a photon in the direction of its electric dipole. This is a well known physical 
property of the dipole radiation approximation in classical electrodynamics [51]. 

Second, the scattering cross section is very small except when the frequency of the incident 
photon coi is close to one of the atomic frequencies co n i ■ The atomic frequencies may correspond to 
energy levels in atoms, for instance, but they may also encompass rotational and vibrational energy 
levels for the case of molecules. 


3.2.3 THE SCATTERING TENSOR 

It is convenient to define the Scattering Tensor as [6]: 


(Cab)fi 


1 y - ' | (Aa)/«(Afo)«( (Afc)/n(Aa)«r 

h „ 1 CO n i — COi CO n f + tWj 


(3.22) 


and use it to rewrite the scattering amplitude in a much simpler form: 


do — U>i(u>i + 0 >if) (C ab )fj (^a) ' (Cfc) 


2 dQ 


(3.23) 


Furthermore, we can decompose the scattering tensor into scalar (0), symmetric ( S ), and 
antisymmetric (A) components: 


Cab — C°8 a b + C S ab T Ca b 


(3.24) 


which are given by: 

C° — -C 

^ ^ '-'aa 

Cab = j ^ ab + Cba) ~ C°8 a b 

Cab = 2 ^ab-Cba) 


(3.25) 
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and summation over repeated indices is implied. Then, the transition amplitudes for the three 
components of the scattering tensor operator are given by: 


(C°) f i 

(C,ab>fi 
(Cab)fi 


tE 


to)ni H“ to)nf 




3 “ (0) ni ~ (Di)((O n f + COi) 

1 (to^ni + to)nf) ((fi'a)fn(fi'b)ni H - (A&)/«(Afl)m) //7 0\ o 

o / v 7 77 i 7 vV ) fi °ab 

2 " 0 to>ni ~ to)i){CO n f + COi) 

2 COi + \ (/^a) fni.l^bini (A/?)//i (Aa)wi 


(3.26) 


E 


iS’hl i ) (OJ/i f 0)j ) 


3 . 2.4 FREQUENCY INVARIANCE 

Let us consider now the scattering of a photon by A identical atoms. The scattering tensor for the 
total system is given by the sum of the individual scattering tensors for each atom. If the atoms are 
close to each other, and all of them have the exact same physical properties, we can expect that the 
scattering tensor for each of the atoms will be the same. The only possible difference in all the terms 
will be due to the initial and final states taken in the transition amplitude. That is: 

(Cabtfr" = (327) 

7= 1 
N 

= e< */ } I I > 

j =i 

Indeed, the wave functions for the atomic states are not the same, but they are defined within 
arbitrary phase factors. Thus, the initial and final states for the j th atom can be written as: 

| vpW j = e i4i | ) (3 28) 

I ) = e^h^f) 

Notice that, in order to avoid confusion with the indices, we use I to denote the imaginary number. 
The total scattering tensor now takes the following form: 

(C ab )f 0,al) = \C ab | Wi) (3.29) 

7=1 

7=1 


As a consequence, to compute the differential cross section we need to take the average over 
the phase factors for each atom separately [6]. However, for the case of Rayleigh scattering, not only 
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the energy is unchanged, but also the phases of the initial and final states are the same: 


and therefore: 


dcr Rayleigh ^ 


oc 


oc 


1>i = 


(Ca) f ( c b) (Cab)^ 


(Total) 


N 


( c a) f *( c bY J2 


]=• 


(3.30) 


(3.31) 


(c c y*(c b y N (C ab )/« 


cx N 2 


(e a ) / *(e fc )‘ (Cab)ft 


On the other hand, for Raman scattering the phase factors for the initial and final states are 
different: 

0/ # 0} 


(3.32) 


and: 


dcr Ra 


oc 


oc 


oc 


oc 


(Ca) f *(c b y (Cab)%° ,a,) 


N 


tea) f *teb)‘ ^ 


1=1 

N 


(3.33) 


7=1 


(,ea) f *(e b y (Cab)fi + tea) 1 * (eb)' (Cab) fi 


/*<WV 


+ 


tea) f *te b y (Cab) fi 


+ crossed terms 


However, all the crossed terms vanish at the time of taking the average over the different phases of 
the atoms [6]. Therefore: 


dcr Raman CX N 


tea)^ teb) (Cab)fi 


which implies that: 


dcr Rayleigh 


N 


(3.34) 


(3.35) 


dcr Raman 

As a consequence, for a large number of atoms, the contributions from Raman scattering are 
much smaller than those that arise from Rayleigh scattering: 


dcr R a yi e jgh dcr Raman for N 1 


(3.36) 
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In the case of radiation scattering, this means that the frequency of the scattered radiation is extremely 
likely to be the same as the one from the incident radiation. Thus, we have obtained the specular 
reflection property P4 solely based on QED arguments. 

3.2.5 QUARTIC FREQUENCY DEPENDENCE 

In the case of Rayleigh scattering, we have that: 


a>i = co f = co (3.37) 

It is easy to verify that the components of the scattering tensor have the following co-independent 
limit values in the low energy regime: 

= 0 (3.38) 

7 ^ 0 
7 ^ 0 


lim (0/7 

co —>-0 

lim S C ab)fi 
lim (C°) fi 

co —>0 


That is, the scattering cross section only depends on the frequency co through the term outside the 
square absolute value of Equation (3.23): 

lim da oc &>,• (co,- + co,j) 3 (3.39) 


and as a consequence: 


lim da oc co 4 

co —>-0 


(3.40) 


This means that, in the low energy limit, the Rayleigh scattering has a quartic frequency depen¬ 
dence [6]. In other words, we have shown that the specular reflection property P5 arises solely from 
QED considerations. 


3.3 MIRROR-PHOTON SCATTERING 

In the preceding section we discussed the scattering of light by a set of atoms. We made no as¬ 
sumptions as to the specific geometry of the set of atoms. However, geometric considerations are 
important to derive the remaining properties of specular reflections. Indeed, the geometry of the 
mirror relative to the incoming photon embodies the angular information required by the law of 
reflections. 

As we will see, the geometry of the scattering system is ultimately responsible for an interfer¬ 
ence pattern at the detection point. To arrive to this conclusion we will model a specular surface as 
a set of scatterers and analyze the interference patterns that they produce. As such, the approach to 
follow is to perform an interferometric analysis similar to the one used in the description of Young’s 
double slit experiment. 



3.3. MIRROR-PHOTON SCATTERING 31 


3.3.1 YOUNG’S DOUBLE SLIT EXPERIMENT 

As a first step, let us consider the standard setup for Young’s double slit experiment (Figure 3.3). 
The slits in the first screen are located at the points rj and r 2 . We are interested in estimating the 
intensity at the measurement point r on the right-hand side screen at time t. 


Y 



Figure 3.3: Diagramatic description of the standard setup for Young’s double slit experiment. 


As seen in Figure 3.3, the electromagnetic field at some point r in the photodetection screen 
has two contributions, one for each slit. The contribution from slit i is the value of the field in r,- 
propagated at the speed of light c across the distance S,- over the time t, \ 

E (+) (r, t ) = E (+) (r ; -, t - t ,-) (3.41) 


where: 



(3.42) 


c c 

From the classical theory of diffraction, it is know that the electromagnetic field at (r, t) 
is given by the linear superposition of both contributions. As a consequence, we have a similar 
expression for the quantum electromagnetic field: 

E (+) (r, t) « Ki E (+) (ri ,t - t{) + K 2 E (+) (r 2 , t - t 2 ) (3.43) 


where the constants K \ and K 2 are pure imaginary numbers that depend on the specific geometry 
of the slits [8, 43]. 

As discussed in the previous chapter, the expectation value of the intensity at the measurement 
point is given by: 

(7(r, t)) = (3.44) 

= ( E (_) (r, t) ■ E (+) (r, t) ) 

which implies, for the case in consideration, that: 

</(r,0> = |X 1 | 2 (E«(ri,t-ti).EW(ri,t-ti)) (3.45) 

+ | K 2 \ 2 ( E (_) (r 2 , t-t 2 )- E (+) (r 2 , t-t 2 )) 

+ 2 Re (k*K 2 ( E ( -)(n, t - ti) ■ E (+) (r 2 , t - t 2 ) )) 
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In terms of the first order correlation functions defined in the previous chapter, the expectation 
value of the intensity is given by: 

(/(r, t)) = |^i| 2 G (1) (n,n;0) (3.46) 

+ \K 2 \ 2 G (1) (r 2 ,r 2 ;0) 

+ 2Re^K*K 2 G m (r l , r2 -,h-t2)) 


If we define: 


(/ ( °(r,0> = \Ki\ 2 G (1) (r,-,r,;0) 

(3.47) 

and: 

dir s G (1) ( ri,r 2 ;r) 

g T r i, r 2 ; r) = = 

^G^(ii,ti-0)G^Hi 2 ,r 2 -,0) 

(3.48) 

where: 


r = t 2 — t\ 

(3.49) 

then we can write: 


(7(r,0) = </ a) (r,0) + </ (2) (r,0> 

(3.50) 


+ 2 yj {/0)(r, r, 0) (g (1) (ri, r 2 ; r)) 

Ifwe perform the polar decomposition of g (1 '(ri, r 2 ; x): 


g m (r l ,t 2 -r) = 


g (l \ ri,r 2 ;r) 


O i(a-o)x) 


where: 


a = a(ri,r 2 ; r) = arg (g (1) (ri,r 2 ; r)) + cox 
and co is the frequency of the field^. Then, we finally obtain: 

</(r,0) = </ (1) (r,t)> + </ (2) (r,0> 


(3.51) 

(3.52) 

(3.53) 


2V</ (1) (r,0></ (2) (r,0} 


^ (1) (ri,r 2 ;r) 


cos (a — cor) 


The cosine of cor is the term that varies rapidly along the surface of the detecting screen, producing 
the interference pattern [92]. Clearly, this term depends on the difference of the distance traveled 
by the light from each of the slits to the detection point: 


co 

cox = — (Sn — Si) 
c 

"’This specific decomposition is taken so that a becomes a slowly varying function of its variables [92]. 


(3.54) 
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Therefore, we can understand S; as a form of interferometric distance. 

So far we have derived the observed interference pattern as a consequence of the variations 
in the source fields Ei(ri,fi) and E 2 (r 2 , t 2 ) at the detection point (r ,t). However, Feynman’s 
path integral formalism provides an alternative interpretation of the results. Within this theoretical 
framework, it is reasoned that the photon could have taken two paths to arrive to the detector point 
(r, /): it either went through (ri, ?i) or through (H, ^ 2 )- The path uncertainty that the experiment 
cannot resolve is what ultimately causes the quantum interference in the screen [29]. 

3.3.2 YOUNG’S DOUBLE SCATTERER EXPERIMENT 

Let us now consider Young’s experiment, but we replace the slits by two scatterers as shown in Figure 
3.4. To simplify the discussion, let us assume that each scatterer is a two level atom. We denote by 
bj and e, the ground and excited states of atom i , respectively. We further presume that a photon is 
scattered by these atoms at a time t such that: 

—8<^t<^8 and 0 < 8 <5C 1 (3.55) 



Figure 3.4: Young’s interference experiment with two scattering atoms. 


For this example, the process of photon scattering is extremely simple [92]. The initial system 
is made of a photon and two atoms in the ground state. The photon is absorbed by one of the atoms, 
which is sent from the ground to the excited state. Subsequently, the atom emits the photon and 
goes back to its ground state. The quantum states of the system before, after, and while the photon 
has been absorbed are given by: 


!*>»<-« 

l v k)-^< f <5 




\y)\blb 2 ) 

A=\ 0 )Qeib 2 ) + \bie 2 )) 
V2 

4 = (In) + In)) \b1b2) 

V2 


(3.56) 


where <5 > 0 is a small time interval. 
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We now define the photon operator yj given by the following property: 

I Yi) = Yi\0) 


(3.57) 


so we can write the state of the system as: 


m t >8 = ~^(y 1+ Y2)\0)\bl b 2) 


(3.58) 


For times t > 8, the expectation value of the intensity measured at the point (r, t ) in the screen 
is given by (dropping the “t > 8” subindex): 

(7(r, t)) = <'P|E ( - ) (M) • E (+) (r,/)|T) (3.59) 

= |<0|(yi + y 2 )E^(r, 0 • E (+) (r, t)(yf + y 2 + )|0) 


(0| [E (+) (r,t), y/l + lE^fr.O, V 2 1 |0> 


^ 1) (Ar l ,t) + *^(Ar 2 ,t) 


»(2 h 


^ I} (A n,t) 


^ 2) (A r 2 ,t) 


■2Re{^\A n j)^f\Ar 2 J)) 


where we have used the photon wave function tyy as defined in the previous chapter, and also: 



(0 [i (+) (r, t), y/ 

0) = (0|(£W(r,Oy i t -y/E (+) (r,o)|0) 

(3.60) 



= (0|E (+) (r,t)y/|0) 

= <0| J E ,(+) (r, t)\Yi) 

= *«(A n,t) 


and: 

'Fy ) (Ar 1 -, t) 

= - A n /c) e -(/-+r/ 2 )( f -A n /c) 

An 

(3.61) 

with: 


> 

II 

i-i 

i 

li 

Co 

(3.62) 

the distance from the i th atom to the point r in the detecting screen. 



As before, the cross term can be rightfully identified as the rapidly oscillating term that causes 
the interference phenomena: 


27?e('lA 1) (Ari, f)'P* (2) (Ar 2 , f)) oc cos (“(Si - S 2 )) (3.63) 

Therefore, the photon interacts with both scattering atoms in a quantum superposition. The 
experiment cannot resolve the atom that scattered the photon, and such uncertainty produces the 
interference pattern measured in the screen. This results in the exact same interference pattern as 
the one seen in the double slit experiment. 
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3.3.3 YOUNG’S MULTIPLE SCATTERER EXPERIMENT 

We can easily generalize the previous results for the case of a Young’s interference experiment with 
N scattering atoms (Figure 3.5). The expectation value of the intensity at (r, t ) is given by: 


</(r, 0) 


1 

N 


A *;.0 

i =1 


(3.64) 


where as before: 

^\ARi,t) = -^-0(f - ARj/c) e -(^+r/2)(t-ARi/c) (3.65) 

y Ar id 

and Arid is the distance from the i th atom to the point r in the detecting screen, and ARj is the 
total interferometric distance of the experiment: 


A Ri = A r si + Artd 

= \r s - T t \ + \ti - td\ 


(3.66) 


and r s , id, and r; are the positions of the source, the detecting point, and the i th atom, respectively. 



Figure 3.5: Young’s interference experiment with N scattering atoms. 

As before, the photon interacts simultaneously with all the atoms and the scattering is mathe¬ 
matically represented by a quantum superposition. The uncertainty of the specific atom that scatters 
the photon is what produces an interference pattern in the screen. 

Notice that in this case, the photon wave functions require an extra phase factor. Indeed, for 
the case of two scatterers it is assumed that the photon source is symmetrically placed between the 
both of them. That is, it is presumed that the distance from the source to the first scatterer is the 
same as the distance from the source to the second scatterer. In the case of three or more scatterers, 
we cannot make such an assumption unless all of them are located in a circumference around the 
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source. The difference of distance from the source to each of the scatterers also contributes to the 
interference pattern. Therefore, the exponential and the step function require the extra phase factor 
to compensate for the right geometry of the scatterers. On the other hand, the amplitude term 
indicates the isotropic decay of the amplitude from its source (the scatterer) to the screen. As such, 
this term does not need to be corrected by the geometry of the scatterers 6 . 

3.3.4 SPECULAR REFLECTION AS A SCATTERING PROCESS 

The preceding discussion suggests that we can describe specular reflections as a scattering process. 
As depicted in Figure 3.6, a light source emits a photon y which is subsequently scattered by the N 
atoms that make the mirror’s surface. Therefore, the detector in the right upper corner will measure 
an interference pattern similar to the one obtained by Young’s experiment with multiple scattering 
atoms. 



Figure 3.6: Specular reflection as a scattering process: a single photon is simultaneously scattered by all 
the atoms forming a quantum superposition (the shape and color of the radiated waves are for illustration 
purposes only, no physical meaning should be given to them). 

In other words, the incoming photon simultaneously interacts with all the atoms in the surface of 
the mirror. This effect is mathematically described by a quantum superposition similar to the one 
used in the analysis of Young’s interference experiment with multiple scattering atoms. Indeed, the 
photon is scattered by the atoms and the intensity of the photonic electromagnetic quantum field is 
measured by the detector. As in the previous examples, there is uncertainty in the actual atom that 
scattered the photon. Such an uncertainty ultimately leads to the interference pattern measured by 
the detector. 

The mathematical description of specular reflection of a photon is the same as the one used in 
the previous subsection. We are only required to understand r a as the position of the photodetector. 

Ahe exact same argument holds true for the case of slits instead of scatterers. 
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As a consequence, for all means and purposes, the description of specular reflections boils 
down to a Young-type experiment with a large number of scattering centers and the photodetection 
screen is replaced by a photodetector on the same side as the light source. Thus, thinking about 
photons interacting with the atoms in a mirror as a sophisticated game of billiards is completely 
wrong. 

Therefore, we have produced a consistent theoretical description of specular reflections solely 
based on fundamental QED scattering processes. That is, we have described specular reflections 
solely in terms of absorption and emission of photons by the mirror’s constituent atoms. 

3.3.5 THE CLASSICAL PATH OF THE PHOTON 

The analogy with Young’s scattering experiments has shown that a single photon is simultaneously 
scattered by all the atoms in a mirror. The scattered photon is described by a quantum superposition. 
The elements in the superposition correspond to the scattering events produced by each of the atoms. 
Uncertainty in the actual scattering event leads to quantum interference. Thus, the detector is really 
measuring an interference pattern similar to the one found in Young’s experiments. In addition, in 
the previous chapter we proved that it is mathematically impossible to localize a photon. It was also 
mentioned that within QED and QFT it is not even possible to consider anything but the initial 
and final states of the system. But nevertheless, the law of reflections seems to require a very specific 
path for the light to travel, which appears to contradict our QED arguments. 

We can reconcile both views by defining the Classical Path of the photom the classical path of 
a photon traveling from a source to a detector is the imaginary path that results from demanding 
maximum constructive interference at the detection point. 

Let us now calculate the classical path of the photon for a specular reflection experiment. 
Imagine the experimental setup shown in Figure 3.7. A single photon light source and a detector 
are located a distance a apart and a distance b over a mirror. Now consider only two of the many 
possible scattering events for a photon reflected off a mirror: the photon is scattered by the atom at 
x — 8 and the photon is scattered by the atom at x + <5. The straight lines in Figure 3.7 represent 
two of the many possible interferometric distances: they do not represent the paths of the photon. 
Indeed, as we discussed before, photon coordinates and photon paths are concepts that do not have 
meaning in the traditional context of QED. 

As discussed in the previous sections, the intensity measured by the detector will have a rapidly 
varying interference term proportional to the cosine of the difference of the interferometric distances. 
That is: 

(Ifd)) oc cos Al ^ A2 ^ (3.67) 

where is the location of the detector and Ai and A 2 are given by: 


Al 

Ai 


h+h 

h +14 


(3.68) 
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Figure 3.7: Two possible scattering events for a photon reflected off a mirror: the photon is scattered 
by the atom at x — 8 and the photon is scattered by the atom at x + 8 . The straight lines represent 
interferometric distances and not photon paths. 


and: 


l\ = b 2 + (x + S ) 2 (3.69) 

l\ = (a — x — 8) 2 + b 2 
l 2 = (x - 8) 2 + b 2 
I 4 = (a — x + 8) 2 + b 2 


Let us now assume that 8 is very small, of the order of the inter-atomic distance. Then: 

IA i - A 2 | » 0(8) =» (I(r d )) cx 1 - (") 2 r(x) (3.70) 

where we have defined: 

tWE(Ai(r)-A 2 (t)) 2 (3.71) 

We need to determine the value of x which produces maximum constructive interference at 
t d . This happens when the variation of r (x) is minimal, as this allows consecutive beams to increase 
the interference amplitude more rapidly. In other words, we need to find the point in the mirror for 
which \dr/dx\ is minimal. 

Figure 3.8 shows the behavior of t(x) and dr/dx when a = b = l,and<5 = 10” 4 .The region 
for which \dr/dx\ is minimal corresponds to the region around x Rs 0.5: 


dr 

dx 


1 

2 


(3.72) 
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Figure 3.8: Maximum constructive interference in two photon scattering events occurs when \dx/dx\ & 
0 =>• x & 0.5. 

Furthermore, the geometry of the experiment shown in Figure 3.7 implies that x & 0.5 is the exact 
same region of the mirror where: 

6 j « e r (3.73) 

Therefore, the classical path of the photon corresponds to an imaginary path where the 
incident and outgoing photons make equal angles with respect to the normal to the mirror. We 
have successfully recovered the law of reflection, specular reflections property PI, exclusively using 
quantum electrodynamical arguments involving the absorption and emission of photons'. 

While the previous example used a one-dimensional mirror, one could introduce another 
distance variable y orthogonal to x to represent the full two-dimensional extension of a flat mirror. 
In this case, the minimum of | Vr(jc, y)| is found in (x, y ) = (0.5, 0.0), which implies that the 
incident classical path of the photon i, the outgoing classical path of the photon r, and the normal 
to the mirror n are all contained in the same plane n. This is specular reflections property P2. 

Finally, it can be observed that the minimal value of \dr/dx\ is reached on a minimum of 
x(x). In other words, the region x & 0.5 not only marks the region where \dx/dx\ £» 0, but it is also 
the region where the distance emitter-mirror-detector is minimal. Therefore, the classical path of 
the photon takes the path of least time. Therefore, Fermat’s principle, specular reflection property 
P3, holds true. 


^ An alternative derivation requires the analysis of all the possible paths taken by the photon [30]. In the context of a path integral 
formulation of quantum mechanics [29], the interferometric distances can be safely interpreted as actual photon paths. Thus, 
Figure 7 would be showing two of the many possible paths that could have been taken by the photon. As discussed in Section 
3.3.1, uncertainty about the traveled path leads to interference. Therefore, the fact that we do not know what atom scattered the 
photon leads to a path uncertainty and to an interference pattern in the detector. 
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3.4 FURTHER DEVELOPMENTS 

In this final section we will address a few natural questions regarding the quantum nature of specular 
reflections. In particular, we will discuss why some materials are better reflectors than others, how 
big or how small a mirror can be, and what is the relevance of the inter-atomic distances and the 
number of incident photons in the specular reflection of light. 


3.4.1 GOOD REFLECTORS 

It is a well known fact from our daily lives that most materials will reflect light to a greater or to 
a lesser degree. However, metallic surfaces are often used to manufacture “good mirrors”. Indeed, 
high quality mirrors are made by applying a reflective coating (silver), to a suitable substrate (glass). 
However, the substrate is only used to protect the reflective metallic coating from corrosion or any 
other type of physical damage. The unique physical characteristic that makes a metal an efficient 
reflecting material is that it has a good conduction band [8]. That is, metals have a partly empty 
electron energy band overlapping the valence band [3]. 

Let us recall that to derive the expressions for the Young’s Scattering Experiments, we have 
assumed the scatterer to be a 2-level atom. Therefore, so far we have limited our discussion to 
reflecting surfaces made of atoms with only two energy levels. That is, these mirrors will only reflect 
light with a frequency near to the atom’s resonant frequency: 


COi COab 


Eg - E h 
h 


(3.74) 


Therefore, so far we have described the quantum behavior of a good reflecting surface for light in 
the close neighborhood of co a b- 

However, the scattering cross section for a photon-atom interaction given in Equation (3.20) 
holds for atoms with an infinite, but discrete number of energy levels {n = 1,2...). We could extend 
these formulas to derive a Weisskopf-Wigner model for these types of atoms. The result would be 
a quantum description of a mirror that is a good reflector for light in the near neighborhood of the 
resonant frequencies a >„. As a next step one could extend this model to the continuum of energies to 
represent the conductivity band of metals. Therefore, the conductivity band of a metal is ultimately 
responsible for the material to be a good reflector over a wide range of frequencies. 

Furthermore, let us recall that the penetration distance d of an electromagnetic wave of 
frequency co incident on a metallic object is given by: 


y/Slt/J.CT a> 


(3.75) 


where o is the conductivity and the magnetic permeability of the material [8]. On the other hand, 
the reflectance of the material is: 


n « i - 2 


CO 


2 n a 


(3.76) 








in the large cr regime [8]. Then, for a perfect conductor: 
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lim d = 0 (3.77) 

G —> OO 

lim 1Z = 1 

G —> OO 

That is, a strong absorber is also a strong reflector [8]. In this situation, the incident photon is very 
likely to be absorbed and re-emitted by the atoms near the surface of the mirror. Very few atoms 
interact with the photon, even though they are strong absorbers. As a consequence, little energy is 
dissipated by the atoms and most of the energy reappears in the reflected light. Thus, for metallic 
surfaces which are good conductors, there is no need to look at the interferometric pattern on the 
“other side” of the mirror. 

In addition, it is also important for a good mirror to be perfectly flat. Otherwise, the light will 
be reflected in many directions, creating the effect of diffuse light. In this case, the law of reflections 
is only valid on an infinitesimally small portion of the surface. 

3.4.2 THE SIZE OF THE MIRROR 

Let us consider the same setup as the one shown in Figure 3.7. What is the minimum size of the 
mirror so we can obtain the law of reflections? If the mirror is much smaller than A,-, the wavelength 
of the incident photon, then a large portion of the incident light will be diffracted. The equations that 
describe the quantum nature of diffraction are exactly the same as the ones used in our discussion 
for reflections. However, in this situation we will not obtain the law of reflections because there is 
not an unique observation point characterized by maximum constructive interference [61]. 

A similar question comes to mind: how big can the mirror be to reflect a single photon? 
As previously discussed, the photon may interact with all the atoms that form an object forming 
a quantum superposition. The size of this range will depend on the optical aperture of the photon 
source. Thus, in principle, there is no upper bound on the size of the mirror. 

3.4.3 INTER ATOMIC DISTANCES 

In the previous subsection we considered the overall size of the mirror. Here we will consider the 
inter-atomic distances between the atoms that make the reflective material. Let us consider first the 
case where the inter-atomic distance is larger than A,-, the wavelength of the incident photon. In this 
situation, the light will be diffracted by each atom and the constructive interference on a privileged 
direction will be greatly reduced [61]. Such a phenomena is known as a diffraction grating and can 
be observed when looking at the surface of a CD or DVD from a low angle (the spiral grooves that 
encode the information form a diffraction grating) [43]. 

On the other hand, there is no lower limit to the size of the inter-atomic distances, and these 
can be as small as necessary (up to the limits imposed by the physical model of QED). 
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3.4.4 NUMBER OF PHOTONS 

In the previous section we described how a single photon gets reflected off a mirror due to scattering 
processes. The law of reflections is recovered by considering the classical path of the photon, which 
is the imaginary line where the constructive interference is maximal. Of course, the interference 
pattern is not made of a discrete set of values, and the intensity has a natural width that goes from 
some minimal to maximal values. 

In the case of two photons, the same ideas apply. However, it can be shown that in this case 
the width of the intensity is narrower [61]. As we continue to increase the number of photons, the 
intensity band becomes narrower and narrower, until it becomes very close to a thin straight line 8 . 
This is not surprising: as the number of photons increases, we approach the classical limit. In the 
limit of an infinitely large number of photons, the photons form well-defined rays of light that are 
reflected on a single direction. That is, the intensity of the field outside the established classical path 
of the photon goes to zero very rapidly. 

3.5 SUMMARY 

We have shown that all five physical properties that describe the specular reflection of light by a 
mirror can be deduced using only quantum electrodynamics and interferometric arguments. That 
is, the processes of emission and absorption of photons by atoms ultimately leads to the law of 
reflection, coplanarity, the Fermat Principle, frequency invariance between incident and outgoing 
radiation, and a quartic frequency dependence of the scattering cross section. 

And no less important, the resulting equations are useful for some applications in quantum 
optics. For instance, as we will see in the following chapters, we can use these results towards the 
development of strategies and procedures to define and estimate quantum radar cross sections [61]. 


8 We will discuss this behavior in Chapter 6. 
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CHAPTER 4 

Classical Radar Theory 


In this chapter we will discuss the theory behind the operation of a classical radar. We will derive 
the radar equation and discuss at some length the radar cross sections that characterize spherical and 
flat rectangular targets. As this is a classical system, we have to assume that the number of photons 
is very large, and as a consequence, classical radar can be perfectly described by Maxwell’s equations 
without recurring to quantum theory. 


4.1 THE RADAR EQUATION 

Figure 4.1 describes the conceptual operation of a monostatic radar system 1 . A transmitter sends an 
electromagnetic pulse towards an airborne target located at a distance R. The power and gain of the 
transmitter are denoted by P t and G t , respectively. The target partially reflects the electromagnetic 
pulse back to the receiver. The power and aperture of the receiver are P r and A r , respectively. The 
pulse takes At from the moment it was transmitted to the moment its reflection is received. 

From the knowledge of At it is possible to estimate the range to the target as: 


R = 



(4.1) 


where c is the speed of light. However, this equation does not express how weak or how strong the 
measured signal will be. 

The strength of the received radar signal is determined through the radar equation [31, 59, 
95, 96]. It can be easily derived through an energetic analysis: if the outgoing power at the radar is 
P r G t , then the incident power density at the target is: 


w, m = P,0, x W; 


(4.2) 


where F € [0, 1] is a form factor that measures the transparency of the space between the radar and 
the target. If this space is perfectly transparent to the radar pulse then F = 1, and if it is perfectly 
opaque then F = 0. For atmospheric attenuation due to fog and clouds, for instance, we can write 
F as: 

F = e~ XcR/2 (4.3) 

where Xc is the attenuation coefficient due to absorption and scattering. 

■*4n a monostatic radar, the transmitter and the receiver are located in the same place. In a bistatic radar, the transmitter and the 
receiver are in different locations. In what follows, we will concentrate our discussion on monostatic radars. 
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Figure 4.1: Basic operation of a monostatic classical radar. An electromagnetic pulse is sent by the radar 
transmitter towards the target. The reflected pulse is measured by the radar receiver. 


The 1 /An R 2 term is the well known decrease of transmitted power of a spherical wave at a 
distance R from the source. We assume that the transmitter is isotropic and radiates uniformly in all 
directions. Thus, the radiated density in a sphere of radius R is given by the radiated power divided 
by the surface area (AnR 2 ). 

Now let us assume that only a portion of this energy is scattered. Then, the scattered power 
at the target is: 

P S @T = Wj@r x a (4.4) 

where a is called the radar cross section and will be discussed in detail in the following sections. Going 
back to our analysis, the scattered power density at the receiver is: 


F 2 

W s @ r = P S @T X 2 

And finally, the scattered power at the receiver is: 


P,G t aF 4 
( An) 2 R 4 


(4.5) 


P r — Ps@r — W s @rA r 

As a consequence, the radar equation is given by: 

P,G t A r crF A 
( An) 2 R 4 

and it is the most fundamental equation in classical radar theory. 


(4-6) 


(4.7) 


4.2 MAXIMUM DETECTION RANGE 

The radar detector is sensitive to thermal noise [53, 95, 96]. The total amount of thermal power Nt 
can be written as: 


Nt = k.BT e B 


( 4 . 8 ) 
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where k B is Boltzman constant, T e is the “effective temperature”, and B is the bandwidth of the 
receiver. 

The signal-to-noise ratio (SNR) is given by: 


SNR = — 

N t 

P,G,A r oF 4 


(4-9) 


(4tt ) 2 R 4 k B T e B 

Now, the aperture area A,- of the receiver can be written in terms of a “receiver gain” G as: 

An A r 


G r = 


and the SNR equation looks like: 


SNR = 


A 2 

P t G,G r X 2 oF 4 


(4.10) 


(4.11) 


(4jr) 3 R 4 k B T e B 

Clearly, there is a minimum value of SNR that allows the detection of a signal. This imposes 
a limit for the detection range of the radar system: 

i 2^ ;r4 \ !/ 4 

(4.12) 


R„ 


/ P,G t G r X 2 (j F‘ 

\(4n) 3 k B T e B SNR m i n 


where in practice, SNR m i n is in the order of 10 to 20 decibels [53]. 


4.3 RADARJAMMING 

Electronic countermeasures are often used in the modern battlefield to defeat the radar systems of 
the enemy [53, 80, 103]. There are passive and active countermeasures. The first kind involves the 
deployment of artificial clutter or chaff, which has a large return that hides the echo produced by the 
weapons systems. The second kind is jamming, which is an active countermeasure. This technique 
involves the active injection of noise to the enemy’s radar. 

The effectiveness of the jamming system can be described mathematically through the jam¬ 
ming equation. Let us assume the radar, target, and jammer system shown in Figure 4.2. The target 
is characterized by its radar cross section a. The jammer emits electromagnetic waves with power 
Pj and gain G j. The radar transmitter has power P, and gain Gt . Finally, the radar receiver has a 
gain G r that in general depends on the angular variables 9 and 0. 

Using similar expressions to the ones used to derive the radar equation, we have that the radar 
receives a jamming signal of power Py 2> given by: 

P r (J) = Wj@ r A r 

(PjGj\ x f ^-GAej,^) 
y47r R 2 J y 47 t 


(4.13) 
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Figure 4.2: Jamming as an active radar countermeasure. 

On the other hand, the target signal power is given by the radar equation: 

pit) _ P,G < G r(°t' <t>t) 

r ~ (4 n)*Rf 

We can use these last two expressions to define a signal-to-jam ratio (SJR) as: 


>(0 


SJR = 


>U ) 

r 

P,G, 

PjGj 


(4.14) 


(4.15) 


(51) 

< 

c 


/N 1 I 

\4tt ) 


'G r (e t ,<p t ) 

G r (dj,4>j) 


The burnthrough range happens when SJR^i 1 and indicates when the jamming signal is able to mask 
the target’s return. That is, the target will be hidden in the electric noise injected by the jammer 
when: 

SJR < 1 (4.16) 

Finally, let us note that in general: 


R j 

Rj RS R t — «: 1 


(4.17) 


which is a benefit for the jammer. But on the other hand, for directional radars pointing directly to 
the target: 

G r (9 ,, 0,) » G r (dj ,4>j)=> » 1 (4.18) 

G r \Pj, <pj) 











which gives an advantage to the radar system. 
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4.4 THE RADAR CROSS SECTION 


The radar equation implies that the radar cross section a is defined by: 


Ps@T 

Wj@T 


(4.19) 


As already mentioned, a represents the fraction of the incident power density that is reflected back 
to the receiver. As such, a is a good measure of the “radar visibility” of the target. Indeed, the larger 
the value of o for a given target, the easier it is to detect it. Similarly, a stealth vehicle is the one with 
a low a. 

An equivalent expression for a is given by: 


4nR 2 P S @T 

= - s— x - 

F 2 P,G, 


(4.20) 


This expression makes clear the fact that a has units of area. A common interpretation of a is the 
hypothetical area of a flat object required to reflect the same fraction of the incident power as the 
target in question. Note, however, that a does not correspond to the projected area of the target, but 
depends on other factors as well. 

A far more formal definition of the classical radar cross section is given by: 


<7 


lim 4nR 2 

R—>oo 




(4.21) 


where E i@T is the incident electric field at the target , while Ei@ r is the scattered electric field as 
measured at the receiver [53, 56]. 

Some important properties of a include: 


Strong Dependencies: It can be proved that a strongly depends on the properties of the target. 
That is, it depends on the target’s geometry (absolute and relative size, shape, and orientation), 
as well as its composition (material properties). 


Weak Dependencies: It can be proved that a is approximately independent on the properties of 
the radar system. That is, it depends very weakly on the strength and the position of the radar 
system. In particular, a is independent of R, the range to the target 2 . 


As a consequence, it is clear that a is a property that (approximately) characterizes a specific 
target, and not the radar system and/or its interaction with the target. Thus, a emerges as a critical 

Tvcn though there is an explicit R~ term in the definition of cr, |E 5 @ r |“ depends implicitly on 1 /R~■ Indeed, |E 5 @ r |^ is the 
electromagnetic energy as measured at the receiver, which has decayed by the isotropic F~ JAtt R^ factor from its value at the 
target. 
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concept for the design of stealth vehicles [53, 56, 70]. In this case, the simulation of a for the 
proposed design of a vehicle will provide a good estimate of its “radar invisibility”. In addition, cr is 
also important to characterize the operational performance and capabilities of radar systems. That 
is, given a radar system, what is the minimum a of a target that it can detect. 

4.5 SCATTERING REGIMES 

Depending on the value of the ratio of the characteristic size of the target a and the wavelength of 
the radar A, three scattering regimes characterize the operation of radar [53, 56, 95, 96]. 

Rayleigh Regime (a <$C A): This is the low-frequency case, and the incident electromagnetic wave 
shows little phase variation over the physical extent of the target. That is, at each instant 
of time, each portion of the target is affected by nearly the same electromagnetic field. In a 
sense, at each moment of time the electromagnetic field is approximately constant over the 
entire physical extension of the target. Then, the incident electromagnetic wave induces dipole 
moments which only depend on the size and orientation of the target. 

Resonant Regime ( a ~ A): In this case, the phase of the incident electromagnetic field changes 
along the length of the target. This means that some of the energy of the electromagnetic field 
is “attached” to the target’s surface, creating surface waves that travel over the spatial extent of 
the target. 

Optical Regime ( a 'Ji> A): In the high-frequency case, collective interactions are minimal and the 
body can be considered made of a collection of independent scattering centers. The total 
electromagnetic scattered field is the superposition of all the individual scattered fields. As a 
consequence, the geometry of the target plays a dominant role on the structure of the scattered 
fields. 


4.6 THE RADAR X-BAND 

The radar region includes the range of frequencies between 3 MHz and 300 GHz. In terms of A, 
this region goes from 1 mm to about 100 m. The applications of radar of most interest, at least from 
a military perspective, are towards the detection and tracking of weapon systems and platforms such 
as missiles and aircrafts. These targets are usually large, ranging from one to a few dozen meters. 
The radar frequencies typically used for missile guidance as well as weather, aircraft, and ground 
surveillance belong to the X-band and range from 8-12 GHz (wavelengths between 2.5 and 3.75 
cm) [95, 96]. For example, airport radar is usually found at 10.525 GHz. This range of frequencies 
is usually known as the Radar X-band because of its classified origins during World War II [85]. 

Then, for applications of interest, the characteristic size of the target a is much larger than 
the wavelength of a radar Xx operating inside the X-band: 


a^> Xx 


(4.22) 
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and therefore, these radar systems operate in the optical regime. 


4.7 SCATTERING MECHANISMS 

Targets of interest, such as missiles or aircraft, are usually complex geometrical objects with inhomo¬ 
geneous electric properties. Therefore, the specular reflection of an electromagnetic wave is only one 
of many possible scattering mechanisms that contribute to the accurate computation of the radar 
cross section. 

Figure 4.3 shows some of the scattering mechanisms that play a role in the determination of o 
for a hypothetical air vehicle [53, 56]. These include specular reflections, creeping and surface waves, 
as well as edge and corner diffraction. The complexity of the target also involves returns due to the 
bouncing of the electromagnetic waves inside cavities or around the shape of the vehicle (interaction 
echo). In addition, there are contributions produced by geometric discontinuities, such as changes in 
the curvature of the surfaces, as well as seam and surface discontinuities due to the physical ensemble 
of the airframe. 



Figure 4.3: Principal scattering mechanisms for a hypothetical air vehicle [53, 56]. 

The most important contributions to the radar cross section for a smooth, continuous, and 
convex target are: 

Specular Reflections: This is the mirror-like reflection of the radar pulse. This contribution follows 
the “Law of Reflections”, that is, the incident and outgoing radiation form the same angle with 
respect to the normal of the surface. 
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End-Region Returns: This contribution is due to the finite size of targets, which often causes 
current discontinuities in the end regions of the body. The result is a sidelobe structure that 
contributes to the scattering away from the specular direction. 

Absorption Effects: Depending on its composition and constitutive properties, the target may be 
able to absorb and dissipate a portion of the incident electromagnetic energy. 

Diffraction Effects: Mostly due to tips, edges, and other discontinuities in the target’s body. 

In what follows, we will consistently ignore absorption and diffraction effects from all of our 
discussions, restricting ourselves to specular and end region returns. 


4.8 SPECULAR AND END-REGION RETURNS 


When an electromagnetic wave impacts a target, electric currents are formed on the surface of the 
body. In turn, these currents radiate electromagnetic energy that contributes to the scattered field. 
However, in practice, the currents that travel the target are unknown and are difficult to estimate 
due to the complexity of the body’s geometry and composition. 

The physical optics approximation can be used to simplify the computation of radar cross 
sections [53,56]. This approximation provides good results as long as the interaction takes place in the 
optical regime; that is, the target is an electrically large body. In addition, this approximation breaks 
down as we increase the viewing angle away from the specular direction. As shown in Figure 4.4, 
using geometric optics we can assume that the current on the shadowed side of the target is zero, while 
it is nonzero in the illuminated portion. We have also assumed an incident plane electromagnetic 
wave. 

The orthogonal area projection A± is the area of the shadow created by the target in a plane 
perpendicular to the wave vector of the incident radiation. 

The current in the illuminated region of the target can be approximated as: 

J 2n x H, (4.23) 


where H, is the incident magnetic field and n is a surface normal [53]. 

The radiated field can be calculated using this current and the Stratton-Chu Radiation Integrals: 

E e (x,y,z ) = f f ■ 9 e M i dV r (4.24) 

£>(■*,: y,z) = J J J J • </> e lkr ' r dv' 

where ?? is the field’s impedance [53, 56]. Assuming aTEM (Transverse-Electro-Magnetic) wave, 
the electric and magnetic fields are related by: 


H(*,y,z) « 


r x E(x, y, z) 
9 


(4-25) 


where it is understood that in the far-zone the radial component of the field is approximately zero. 
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Figure 4.4: Interaction of a plane electromagnetic wave with a target in the physical optics approximation: 
using geometric optics we determine that the current in the shadowed area is zero, and non-zero in the 
illuminated portion of the body [53]. 


4.8.1 SPHERICAL TARGET 


As a first example, let us consider a monostatic radar illuminating a spherical target of radius a as 
shown in Figure 4.5. The current created by an incident wave of amplitude //q and momentum k is: 


where: 


h = 2 H 0 e 


v = 
w = 
A = 


,kA (yw — zt>) 

(4.26) 

sin 0 sin </> 

(4.27) 
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a( 1 — w) 



using the angles from standard spherical coordinates 6 and (j) [53]. 

Using the radiation integrals we can express the scattered electric field as: 
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Figure 4.5: Geometry of a spherical target and an incident electromagnetic wave. 


and therefore, the radar cross section is given by: 


a 


= Tta 


2 



sin (ka)\ 
ka ) 


(4.29) 


which clearly manifests a spherical symmetry (functionally independent of the angular variables 0 
and 0) [53]. 



Figure 4.6: Radar cross section a for a sphere of radius a and wave-vector k. 
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The dependence of a on ka is shown in Figure 4.6. Two different regimes are observed. The 
low-frequency Rayleigh Regime for which: 


2 (kdy 4 

ka < 3 C 1 =>■ a na — 7 ^- oc co 


36 


and the high-frequency Optical Regi, 


ime: 


ka 1 =£• a & na 1 = A_l 


(4.30) 

(4.31) 


Notice that in this regime, a is equal to the projected area of the target Aj_. This is a specific result 
for spherical bodies, and in general, for other geometric objects, this relation will not be true. 

However, it is important to state that the Rayleigh regime solution is only an approximation, 
as we have used physical optics methods to derive cr. An exact solution involves a more sophisticated 
technique known as Mie Theory [ 8 ]. 

Furthermore, in the limit of the optical regime: 

lim a = na 1 = A± (4.32) 

which is independent of A. 


4.8.2 RECTANGULAR TARGET 

Let us now consider a flat rectangular plate of size A = a x b as the one shown in Figure 4.7. If 
we only consider specular returns, then the radar cross section is zero everywhere except for 0 = 0 . 
Therefore, avoiding specular reflections for the design of stealth vehicles simply requires that all 
surfaces should never be perpendicular to the propagation vector of the incoming radar pulse [70]. 
This is one of the reasons why the first stealth aircraft, the Lockheed Martin F117, resembles a 
diamond made of flat plates [ 86 ]. This design avoids specular returns in case an enemy ground radar 
is positioned in the forward direction. 

Unfortunately, such a simple behavior of the radar cross section is correct only if the plate is 
of infinite size. Indeed, the finite size of the target introduces end-region returns. Let us consider 
an incident electromagnetic wave that generates a current in the surface of the target, as shown in 
Figure 4.7. Because the plate is a finite body, the current is discontinuous: the current is non-zero on 
the illuminated surface of the target, but it is zero on the dark side of the plate. Such a discontinuity 
is responsible for a contribution to the electromagnetic field that is radiated in all directions. 

We can use the physical optics approximation and the Stratton-Chu radiation integrals to 
quantify the effect of the end-region returns on the radar cross section of the rectangular target [53]. 
First we assume an incident 6-polarized plane electromagnetic wave arriving from an angle (4,0): 

Ef = E oe 0 e~ ih t (4.33) 

If the plate is uniform and isotropic, the surface current j. v will be given by: 

„ g-'k-r 

jv -2 z x E oe 0- 

1 7 


(4.34) 
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Figure 4.7: Sensor-target geometry for a rectangular 2D target. 



Figure 4.8: End-Region returns contribute to the oq- These returns are attributed to the current dis¬ 
continuities due to the finite size of the target. 


The far-zone scattered field will be given by: 


Ee(r, 0, 0) = e~ jkr 
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(4.35) 
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where: 


it = sin 6 cos 0 
v = sin 6 sin 0 


and the radar cross section a is calculated to be: 


cr ( 9 , 0) = 


Ait A 1 


A 2 


■ cos 0 


sin kau 
kau 


2 2 sin kbv\^ 
kbv ) 



(4.36) 


(4.37) 


Figure 4.9: Simulation of ctc vs. 0 (0 = 0) for a 2D flat rectangular plate. The highest value of ere is 
reached when the target is oriented in the specular direction (0 = 0). The sidelobe structure is entirely 
due to end-region returns. Diffraction and absorption effects have been ignored. 

Figure 4.9 shows the plot of cr/A 2 for 0 = 0 and two plates of size ( a = 10A, b = 10A) and 
(i a = 5A, b = 5A). In both cases, the highest peak is found at 0 = 0 and corresponds to the return 
due to specular reflection: 

cr(0, 0) = —— oc co 2 (4.38) 

A 

Thus, the radar cross section has a quadratic dependency on the frequency of the incident electro¬ 
magnetic wave. Furthermore, if we recall that we are working in the optical regime: 

eft 

cr(0, 0) oc — a 2 (4.39) 

A^ 

that is, in the specular direction, the radar cross section is much larger than the physical area of the 
target. 

In addition, Figure 4.9 also makes evident a crucial property of end-region returns: they 
produce a sidelobe structure around the specular reflection region [53,56]. This is a direct consequence 











































56 4. CLASSICAL RADAR THEORY 


of the oscillatory behavior of the angular variables in the expression for the radar cross section, which 
in turn is due to the discontinuities in the surface currents. 

The radar wavelength and size of the target determine how many lateral sidelobes are present 
in the radar cross section. Let us consider the 4> = 0 case. Then, the radar cross section is proportional 
to: 


n I sin (kau) 

cr(9, 0) oc cos" 9 (- 

kau 


(4.40) 


The radar cross section reaches the value of zero on both sides of each sidelobe. This only happens 
when: 

sin (kau) = 0 (4.41) 


which means that: 


kau = 


2tt a 

-sin 9 = 0,7T, 2jt, ... 

X 


(4.42) 


However, the value of 0 corresponds to the maximum of a. Then, the radar cross section is zero 
when: 

— sin 9 e Z - {0} (4.43) 

X 

Then, for example, if a = 2.5X, we have five zeroes in: 


sin# = 


1 2 3 4 5 
5’ 5’ 5’ 5’ 5 


(4.44) 


and therefore, a has four sidelobes on each side of the maximum. In general, the number of sidelobes 
per side is given by: 

[_2 a/X - 1J (4.45) 

It is important to remark that the sidelobe structure is the most important return after specular 
scattering. Indeed, the sidelobe structure can be used to detect targets even if the specular returns 
contribute in a negligible manner to the received radar signal [53, 56, 70]. Therefore, the sidelobe 
structure is crucial for the design of stealth vehicles (how to reduce it) and radar systems (how to 
exploit it). 


4.8.3 GEOMETRY, FREQUENCY, AND SIDELOBES 

As previously discussed, the sidelobe structure of a can be exploited to enhance detection, even 
if the target avoids specular returns. However, the sidelobe structure can be reduced by smoothly 
tapering to zero the currents over several wavelengths [53]. For example, the sidelobe structure can 
be significantly reduced by changing the conductivity of the plate so that the induced current is 
linearly distributed in such a way that j. v = 0 along the edges. 
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Furthermore, the sidelobe structure depends on the geometry of the target. Let us first consider 
the case of a flat rectangular target. Assume three targets of the same area but different proportions: 

Ai = 2.5 x 4.0 X 2 (4.46) 

A 2 = 5.0 x 2.0 X 2 
A 3 = 10.0x1.0 A 2 

The cross section for these targets is shown in Figure 4.10. From these plots, it is evident that the 
three targets have the same maximum value of a at (9, 0) = (0, 0), but their sidelobe structure is 
different. Therefore, on this specific orientation (0 = 0), target A 1 is more likely to be detected by 
its sidelobe structure than target A 2 . 



Figure 4.10: Radar cross section for a flat 2D plate of equal size, but different proportions 
(Ai:red; A 2 :green; A 3 :blue). 


Another way to reduce the sidelobe structure is by changing the curvature of the surfaces 
that make the target. In double curved surfaces, for example, the currents smoothly taper to zero 
in the shadow regions [56]. The relationship between a and geometry at the specular direction 
(6 = 0 = 0) is illustrated in Table 4.1 [56]. 

It is important to note that the maximum value of the radar cross section has a different 
functional relationship with the frequency depending on the geometry of the target. It is only in the 
case of a doubly curved surface that the radar cross section is independent of the frequency. 

We can use these values to compare the radar cross section of a spherical target: 

<j 0 = no 2 (4.47) 


and a flat square target of similar size: 


47T (2a) 4 
X 2 


(4.48) 
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Table 4.1: Dependence of the radar cross section on target’s geometry [56]. 

Surface 

Example 

Curvature Radii 

a (0,0) 

&> Dependency 

Flat 

Plate 

(oo, oo) 

4jtA 2 /X 2 

2 

(X CO 

Singly Curved 

Cylinder 

(a, oo) 

2nal 2 /A 

(X CO 

Doubly Curved 

Spheroid 

(a\,a 2 ) 

jza\ci2 

(X 1 


in the optical regime: 

2 

— = 64-^ » 1 =>- <t<> » (To (4.49) 

cr 0 A- 

That is, the radar cross section of a sphere is much smaller than the radar cross section of a flat plate of 
similar dimensions. This result is to be expected, as the plate will reflect most of the electromagnetic 
energy back to the radar, while the sphere will scatter a large part of it to other regions. 

Let us consider a specific example. Assume a flat square plate of size a = n'k with A = 0.1 m 
and?; = 100. That is, a = 10 m. Then, the area of the plate is A = a 2 = 100 m 2 . Then, the maximal 
value of its radar cross section at the specular direction is: 



1.3 x 10 7 m 2 


(4.50) 


This means that: 

— ^1.3xl0 5 (4.51) 

A v ' 

which means that the radar cross section at the specular direction is 100,000 times larger than the 

physical area of the plate. 

We can now compare the radar cross section of this target with the one of a sphere of radius 
a with cross section given by: 

(j 0 = 314 m 2 (4.52) 

a z 

As seen in Figure 4.11, the radar cross section of the plate is larger than the radar cross section 
of the sphere for a viewing angle of about 0.31 radians with respect to the specular direction. 

The dependency of the radar cross section on the geometry of the target is one of the reasons 
as to why the new generation of stealth aircrafts does not have the diamond shape of the FI 17. For 
instance, the stealthier B2 Spirit features double curved surfaces across its fuselage [1, 86]. 


4.9 ENERGY CONSERVATION IN THE OPTICAL REGIME 

If we ignore absorption effects, it is reasonable to require that energy is conserved during radar-target 
interactions. In the optical regime it is possible to find a simple equation that neatly describes this 
fact. Although it is not an important result in the context of classical radar theory, this equation will 
prove essential for the simplified description of a quantum radar cross section. 
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Figure 4.11: Radar cross section for a rectangular plate (blue) and a sphere (red) of similar dimensions. 
The former is larger than the second over an angle of about 0.31 radians around the specular direction. 

Let us assume we have the expressions for the scattered fields under a bistatic radar. The 
coordinates of the transmitter are r v and the coordinates of the receiver are r d with respect to the 
target. If the energy is conserved, then the total incident energy on the target has to be scattered 
somewhere in space. That is: 
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IV, (r v ) dS ^ lim 


71(0,0) 


X 


rd ^°° JSr d DT 


W s (r s ,r d ) r j cm d 


(4.53) 


where the first integral is taken over the orthogonal projected area of the target 7j_, and the second 
integral over S rd , a sphere of radius r' that surrounds the target 71 That is, we are integrating over 
all the possible positions of the radar receiver. In terms of the electric fields: 
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'71(0,0) 

If we assume a harmonic field: 

we get: 
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(4.54) 


(4.55) 


(4.56) 


where A± is the area of the orthogonal projection of the target and will depend on the orientation 
of the target. Energy conservation then implies that: 


\E i (r s )\ I A ± (e,<P) 
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(4.57) 
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Using this expression for the conservation of energy we can rewrite the radar cross section for 
a monostatic radar as: 


a = lim 47rr‘ 

r—>oc 


\E S \ 

\Et \ 2 


where we have set: 


: 47r Aj_ lim — 

^°° Is , 


r = L = r d 


\E s {r)V 


\E s (r s ,r d )\ 2 dQ. d 


(4.58) 


(4-59) 


on the numerator. Notice that, in general, the projected area is a function of the angular coordinates 

A L = A L (e,4>). 

This expression for the radar cross section is only valid in the optical regime. Indeed, all the 
incident radiation is absorbed and scattered by a target only if the target’s dimensions are greater 
than the wave length of the incident electromagnetic wave. One way to think about the physics 
behind this fact is to consider the radar as a diffraction limited optical sensor. The angular resolution 
6 is given by the Rayleigh criterion: 

sin 9 rs 1.220— (4.60) 

D 

where D is the diameter of the lens’ aperture and the numerical factor is derived from the position 
of the first dark ring surrounding the central Airy disc [8, 43]. In the optical regime, k and 0 are 
very small, which means that the device can resolve small variations in the geometry of the target. 
On the other hand, if k is large, then the angular resolution is large and may even be larger than the 
target itself. In this case, not all the incident radiation is being scattered, as a portion of it is “passing 
through” the object. 

As an example, let us consider a spherical target Sj of radius a. Then, integrating the above 
expression for a we get: 

[ a JO ~ 477 A_l (4.61) 

JSrDSr 

In addition, because of the spherical symmetry of the target, we also know that: 
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a dQ = g 477 


and therefore: 
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(4.62) 


(4.63) 


which is the value we had already estimated for the spherical target in the optical regime. 

Let us consider now the case of a flat square target of area A = a 2 . From the previous section 
we know that the scattered field is given by: 


\E s (r,e,<j))\ 2 = A 2 


477 “/*" 


E 2 cos 2 9 


sin (kau) 
kau 


sin (kbv) 
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(4.64) 
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and as a consequence, using the traditional expression for the cross section at <p = 0 we obtain: 


47tA 2 0 /sin(A:aM)V 

a(e ,„ ) = _ cos2 ^__) 


(4.65) 


On the other hand, using the expression for the scattered field in the new expression for the cross 


section we have: 

a f sin( kau)\ 

4?rAj_ cos 6 {~l^r) 

cr(6, 0) = —--^-4--- 

I cos2 ° (^) (^) da 


(4.66) 


As we are interested in comparing the accuracy of our approximation, we have added a “tilde” symbol 
to distinguish between both of them. Furthermore, recalling that the maximum value of the cross 
section is reached at 6 = 0, we define the normalized error between both expressions as: 


m = 


om-°m 

<X( 0) 


(4.67) 



Figure 4.12: Relative error using the expression for a based on energy conservation on the optical regime. 
The error diminishes as we get further inside the optical regime. 

Figure 4.11 shows the behavior of e(0) for a plate of size ka = 20n and for a larger plate of 
size ka = 200iz. As the value of ka increases, the error decreases. As expected, the error decreases 
as we penetrate the optical regime (ka 1). 

4.10 RADAR AS AN INFORMATION CHANNEL 

Formally, we can describe radar performance as an information channel. The target information 
produced by the radar system includes [53, 56, 95, 96]: 

Range: The range to the target is obtained through a measurement of the time interval for the echo 
signal At. 
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Speed: The speed of the target can be obtained by analyzing the Doppler shift or changes in the 
target’s range. 

Direction: A directional antenna can provide information about the target’s direction. 

Size: The size of the target has an influence on the magnitude of the returned signal due to the 
radar cross section. 

Shap e: The shape and geometry of the target determine its radar cross section. 

Components: The components of the target, such as moving parts in the case of a propeller, can 
be determined by analyzing the frequency components of the Doppler signal. 

The fidelity of the information depends on the amount of noise in the signal, but also on the 
value of the radar cross section of the target. Therefore, the process of obtaining information from a 
target using radar, can be described as a noisy communication channel. Furthermore, the radar cross 
section could be treated as a form of “geometric noise”. 

4.11 SUMMARY 

The radar equation describes the strength of the electromagnetic signal that arrives to the receiver. 
This equation depends on the radar cross section, which is an important quantity that characterizes 
the visibility of the targets as seen under a radar. Due to the finite size of the targets, a sidelobe 
structure is observed on the angular dependence of the radar cross section. Strongly dependent on 
geometry, these sidelobes enable the detection of targets that avoid specular returns, and therefore, 
they are crucial for the optimal design of stealth vehicles and radar systems. 
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CHAPTER 5 


Quantum Radar Theory 


As discussed in the previous chapter, traditional radar systems employ classical electromagnetic waves 
to detect targets. That is, the traditional radar signal is made of a large number of photons. Quantum 
radar generalizes the concept of radar, but this novel device operates with a relatively small number 
of photons. As such, its theoretical description has to be done using quantum electrodynamics. 

The motivation for developing the concept of quantum radar is the potential use of quantum 
phenomena to increase the sensitivity of the sensing system. As we will see in this chapter, entangled- 
photons quantum radar appears to offer higher sensitivity. Furthermore, it seems that it is much more 
difficult to conduct jamming on a quantum radar system. In addition, quantum radar inherits some 
of the advantages of traditional radar systems. For instance, microwaves photons are able to penetrate 
clouds and fog, while photons in the visible regime are most likely to be absorbed or scattered by 
the environment. 

In this chapter we will present the theory behind quantum radar. We will place emphasis on 
the two major standoff sensing techniques proposed to date: the interferometric quantum radar and 
quantum illumination. Even though it does not operate on the microwave regime that characterizes 
radar, we will briefly discuss quantum LADAR. In addition we will overview the state-of-the-art 
on the technologies required for the construction of quantum radar transmitters and receivers. 

5.1 STANDOFF QUANTUM SENSORS 

In general, we can define a quantum radar as a standoff detection system that uses microwave photons 
and exploits some form of quantum phenomena to enhance its capabilities to detect, identify, and 
resolve targets of interest. 

The targets are assumed to posses a low reflectivity and are far away from the detector. Fur¬ 
thermore, the radar-target system is presumed to be immersed in a noisy and lossy environment. 
That is, the radar signal may be attenuated through absorption or scattering processes, and the overall 
performance of the system is affected by the presence of noise. 

Furthermore, as its name suggests, these devices operate in the microwave regime that charac¬ 
terizes classical radar signals. This is an important feature that inherits the advantages of traditional 
radar pulses such as being able to penetrate clouds and fog. 

5.1.1 THE FUNDAMENTAL LIMITS OF QUANTUM METROLOGY 

The measurement process in quantum mechanics is restricted by Heisenberg’s uncertainty princi¬ 
ple [44, 60]. Suppose we have a quantum system made of N photons and the quantum state has a 
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relative phase 0. Then, in terms of the state phase 0 and the number of photons A, the Heisenberg 
principle leads to the inequality: 

A0 AN > 1 (5.1) 

where A0 and A A represent the phase and photon number fluctuations, respectively [68]. 

Let us also assume we have a monochromatic coherent state of the form 1 : 


i«> = e - |a|2/2 E 

N =0 


„N 


( A !) 1 / 2 


I n) 


(5.2) 


where | n) represents the state of n photons with frequency oo [68, 73, 92]. Then, the probability of 
finding N photons is: 

P{N) = |<A|cx>| 2 (5.3) 

e~W(N) N 

n\ 

which clearly obeys a Poisson distribution law [68]. Furthermore: 



( N) = (a A a) = (oc\ci^ d\a) = a 2 

(5.4) 

and: 

(A 2 ) = (o'|a t aa t a|a) = (A) 2 + (A) 

(5.5) 

And therefore: 

(AN) 2 = (A 2 ) - (A) 2 = (A) 

(5.6) 

which means that: 

AN = ^(N) 

(5.7) 

This is known as shot noise 

In the case of shot 

and it is due to the discrete nature of light in the quantum realm [10,57,68]. 
noise, the Heisenberg inequality takes the form: 


1 

A0 > _ 

vwy 

(5.8) 


This expression for the limit on the sensitivity for a phase measurement is known as the Standard 
Quantum Limit 2 [10, 35, 57]. 

The standard quantum limit is a consequence of the discrete nature of the electromagnetic 
field, as well as the Poissonian statistics of classical light [57, 68]. These properties of quantum light 
translate into vacuum fluctuations that affect the measurement of the amplitude of an electromagnetic 
field [14, 15]. In a sense, these features of the quantum electromagnetic field tend to prevent the 
cooperative behavior of the photons [35]. 

^ Coherent states of light represent quite well the light emitted by a laser. 

2 This limit is also known as the shot-noise limit or coherent state limit. 
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On the other hand, one could argue that the basic principles of quantum theory do not 
impose a specific value on the fluctuations on the number of photons. The standard quantum limit 
is a consequence of the use of non-optimal measurement strategies, and in principle, it could be 
avoided: 

(A 2 N) t* O (<A0 2 ) (5.9) 

That is, the true fundamental limit is the one given by Heisenberg’s principle, which leads to the 

Heisenberg Limit'. 

Ad> > — (5.10) 

(N) 

The Heisenberg limit does not depend on the measurement strategies and it is unavoidable [35, 57]. 

The sensitivity of most of our modern sensors is bounded by the standard quantum limit. 
However, there is reason to believe that it is possible to build sensors that can beat the standard 
quantum limit but are limited by the Heisenberg limit [35, 36]. A regime of variables for which 
the sensitivity of a sensor exceeds the value imposed by the standard quantum limit is called the 
Supersensitivity Regime, and the sensor is said to be Supersensitive. 

In particular, it has been shown that entanglement offers a way to bypass the barrier imposed 
by the standard quantum limit [35, 36]. But nevertheless, it is important to note that entanglement 
is not strictly necessary to achieve the Heisenberg limit. To date some other alternatives have been 
proposed in the literature. These include, for example, the injection of squeezed states [7, 104] or 
the injection of Fock states into the measurement system [46]. 

In any event, most strategies are based on the exact same principle: perform collective mea¬ 
surements after highly correlated states have been injected into the system [35, 36]. Such is the case 
of the interferometric quantum radar that uses highly entangled states to reach the Heisenberg limit. 

Finally, it is important to remark that Heisenberg limit metrology is very challenging to 
accomplish [35, 36]. At these high precision scales, there is a large number of factors that have 
important contributions and cannot be neglected. For instance, it is crucial to estimate the effect of 
thermal noise and other dissipative processes on the signal, as well as any other sources of noise that 
may limit the performance of the detector and the transmitter. 

5.1.2 CLASSIFICATION OF STANDOFF QUANTUM SENSORS 

Standoff quantum sensing architectures can be classified according to the type of quantum phenom¬ 
ena exploited by the system [42]. The three basic categories are the following: 

Type 1 : The quantum sensor transmits un-entangled quantum states of light. 

Type 2: The quantum sensor transmits classical states of light, but uses quantum photo-sensors to 
boost its performance. 

Type 3: The quantum sensor transmits quantum signal states of light that are entangled with quan¬ 
tum ancilla states of light kept at the transmitter. 
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The principal examples of quantum radar systems are the single photon quantum radar and the 
entangled photon quantum radar. These are Type-1 and Type-3 sensors, respectively. On the other 
hand, quantum LADAR is an example of a Type-2 sensor. The major features of these examples 
will be described next. 

5.1.3 SINGLE-PHOTON QUANTUM RADAR 

Single photon quantum radars are Type-1 sensors. These systems operate in a similar fashion as a 
classical radar. As shown in Figure 5.1, a quantum radar transmitter sends a single photon pulse 
towards a target. The target reflects the photon which is subsequently detected by the radar receiver. 

One of the advantages of quantum radar is the fact that some targets appear to look bigger 
when looked with low photon number pulses. That is, as we will see in Chapter 6, the radar cross 
section near the specular region of some targets is larger than its classical counterpart when the 
targets are illuminated with individual photons (entangled or not). 



Figure 5.1: The basic concept of a single photon quantum radar: a single photon is emitted towards a 
target, and subsequently the photon is reflected back to the receiver. 


5.1.4 ENTANGLED-PHOTONS QUANTUM RADAR 

The greatest benefit of quantum radars is obtained with the use of entangled states of light. These 
are Type-3 sensors. As shown in Figure 5.2, an entangled pair of photons is generated. One of these 
photons is sent towards the target and the other kept in the radar system. The outgoing photon is 
reflected by the target and subsequently received by the radar. At this point the correlations embedded 
in the entangled states are exploited to increase detection performance. 
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Examples of these systems are the Interferometric Quantum Radar and Quantum Illumination. 
Both of these are discussed in some detail in the following sections. 



Figure 5.2: An entangled-photons quantum radar system: an entangled pair is produced, one of the 
photons is sent towards the target while the other is kept in the system. The correlation between these 
photons can be exploited to increase the performance of the device. 


5.1.5 QUANTUM LADAR 

Quantum LADAR is a Type-2 standoff quantum sensor that operates in the visible and near-visible 
regimes ’ [22, 23, 41, 58, 79, 93]. As such, quantum LADAR does not count as a form of quantum 
radar technology. Nevertheless, it is important to mention that this is another promising quantum 
sensor which theory is well understood and the development of an operational prototype is currently 
underway. The advantages of LADAR over radar resides on the small wavelengths of the former, 
which allows it to resolve small spatial information [76]. On the other hand, LADAR cannot 
penetrate fog or clouds as radars do, and as a consequence, the operational range of a LADAR is 
usually restricted within 100 km [76]. The use of squeezed vacuum injection to reduce the noise 
due to vacuum fluctuations is one way in which quantum effects can be harnessed to improve the 
performance of LADAR [58]. 


^ LADAR stands for Laser Detection And Ranging, and it is the exact same technology as LIDAR (Light Detection And Ranging). 
However, in the military context, the former is the preferred acronym. 
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5.2 INTERFEROMETRIC QUANTUM RADAR 

The quantum radar systems described in Figure 5.1 and 5.2 resemble an interferometer with a very 
long arm (as big as the range to the target R). This suggests that one way to analyze the performance 
of a quantum radar is to describe it as a quantum interferometer. 


c 



Figure 5.3: Mach-Zender interferometer with input ports A and B, output ports C and D, and a phase 
delay 0 in one of the arms. 

Let us consider the Mach-Zender interferometer made of two mirrors and two beam splitters 
as shown in Figure 5.3. One of the arms of the interferometer introduces a delay phase 0.The value of 
0 can be estimated by measuring the intensity (i.e., number of photons) of the two output beams. In 
agreement with our previous discussion, it can be shown that, after measuring the output for a beam 
with N non-entangled photons, one can estimate the value of 0 with a statistical error proportional 
to 1 /yf~N [20, 66]. That is, this interferometric measurement is bounded by the standard quantum 
limit. 

Similarly, let us assume we want to use an interferometer-like sensor to determine the range 
of a target. The best classical strategy to perform a range estimation is the one currently exploited 
by traditional radar systems. That is, the direct calculation of the average travel time measured for 
N non-entangled photons. It can be shown that in this case the error is given by: 

SR^o( - l —A (5.11) 

Va coVn) 

where Aco is the bandwidth of the radar signal [35]. 

However, as we will see next, entangled photons can be used to beat the standard quantum 
limit and reach the Heisenberg limit. 
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5.2.1 QUANTUM INTERFEROMETRY 

Let us assume we have a NOON state: 


1**00*) = —7= (l-YO) + \0N)) 

V2 


(5.12) 


which is a highly entangled state of two photonic modes in an interferometer [9], These states are 
often proposed for a variety of quantum sensing applications. Each half of the entangled state passes 
through a different arm of the interferometer, so we can write the state as: 


I^noon) = -t= (|AOi|0) 2 + |0)i|Y) 2 ) 
V2 


(5.13) 


where the subindex identifies the arms of the interferometer. Furthermore, we assume that we are 
interested in measuring a phase difference 0 that arises when photons are sent through the second 
arm of the interferometer [20, 66]. The effect of a phase shift 0 on the NOON state is given by: 
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We can write the NOON state in terms of the creation and annihilation operators as: 
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(5.15) 


In order to measure the phase shift 0, the detector has to implement a measurement of the 
following observable: 


A d = |M))(0JV| + \0N){N0\ 

= ((fij^ioxoKfi^ + (4)*|0)(01$)*) 

With this setup, the amount of “noise” is given by: 
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while the phase responsivity is: 
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Then, the phase estimation error is approximately given by: 
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That is, this interferometric phase measurement procedure using highly entangled states can reach 
the Heisenberg limit [20, 66]. 


5.2.2 ATTENUATED QUANTUM INTERFEROMETRY 

Let us assume now the exact same interferometric experiment described above, but the interferometer 
is immersed in an attenuating medium [32, 33] . The expression for the propagation of an attenuated 
NOON state was given in Chapter 2 and is repeated here for convenience: 

1**00*) = (af^lOhlOh (5.20) 

+ —^ e - imr l2/^-Xc2(.co)/2KNL 2 ) |0)i|0) 2 (5.21) 

V2A! v ' 

+ m 


where each component travels across distances L\ and L 2 over two different media characterized 
by refraction indices r]\ and 772 , and attenuation coefficients Xcl and Xc2- Also, |<t>) is a state that 
represents those states scattered outside the NOON basis. 

In the case of the interferometric experiment described above, the amount of noise is given 
by: 

A 2 Ad = - (a^ — la^a^ + <*2 + ( a \ a l) N sin 2 A! (0 — 0o)) (5-22) 

while the phase responsivity is: 
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and we have defined the transmittance of the medium as: 

— 0~XciLi 


and the dispersion phase shift as: 
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Then, the phase estimation error is given by: 
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It is only in the limit of no attenuation: 
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that the Heisenberg limit can be achieved. 
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For comparison purposes it is useful to know the phase estimation error in the case where one uses 
separable states instead of entangled NOON states [32, 33]. The state (without attenuation) is taken 
to be: 

l*s> = -7=(|10> + |01»®* (5.28) 


which can be rewritten as: 


N 


«=^ n («; 

VZ 7=1 

The observable for separable states is: 


W) +4 a) 


3 ) |0)i|0> 2 


(5.29) 


N 

A s = 0 (|01>°' )0) (10| + |10)°' )O) (01|) 


(5.30) 


7=1 


Introducing the attenuation of the states as we did before, we obtain the following variance on the 
phase measurement: 


2 N . 2 

A"As = — (a i — 2a\U2 + 012 ) + N a\_a 2 sin (0 — 0 0 ) 


and a phase responsivity: 


d(A D ) 

dcj) 


= — N Ja\ot 2 sin(0 — 0o) 


As a consequence, the phase estimation error is: 


50 = 


2^1 + 2 ^ - 1 + sin2 (0 - 0 o ) 


*/N I sin(0 - 0 O )| 
which in the no-absorption limit takes the value: 


lim lim 50 = —— 
ai->l a 2 ->-l yiV 


(5.31) 


(5.32) 


(5.33) 


(5.34) 


As expected, separable states cannot be used to beat the standard quantum limit even in the absence 
of attenuation. 


5.2.4 ATMOSPHERIC QUANTUM INTERFEROMETRY 

Typical examples of atmospheric attenuation caused by fog or clouds are given in Table 5.1 for 
different wavelengths in the radar region [95, 96]. In particular, A = 3.2cm is inside the radar X- 
band (2.5-3.75 cm; 8-12 GHz). It can be observed that in this regime, the attenuation coefficient 
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increases as we decrease the wavelength. Therefore, it is convenient to be on the higher-end of the 
radar X-band. 

As a comparison, the average atmospheric attenuation for visible light at 550 nm is of about 
0.87 dB/km, much higher than the attenuation suffered by X-band radar in the cases considered in 
Table 5.1 [91]. This is to be expected, the ability of radar waves to penetrate fog and clouds is the 
reason as to why radar is so important in the military and civilian worlds. 


Table 5.1: Attenuation (dB/Km) caused by fog and clouds 


on different wavelengths in the radar region [95, 96]. 


Visibility (m) 

A = 1.25 cm 

A = 3.2cm 

X = 10cm 1 

30 

1.25 

0.20 

0.02 

90 

0.25 

0.04 

0.004 

300 

0.045 

0.007 

0.001 


In general, the dielectric constant of a medium varies with temperature, and as a consequence, 
the attenuation coefficient varies with the temperature [95, 96]. The values reported in Table 5.1 
for the attenuation coefficient are taken at 0°C. Approximately, the values for 15°C and 25°C can 
be obtained by multiplying the coefficient in the table by 0.6 and 0.4, respectively. That is, in this 
specific regime, the attenuation coefficient decreases as the temperature increases. 

Let us analyze the behavior of the quantum interferometer in the case where A = 3.2 cm. 
We consider two of the scenarios shown in Table 5.1: a visibility of 300 m which is equivalent to 
an attenuation constant of Xc = 0.0016/ km, and a visibility of 30 m equivalent to an attenuation 
constant of Xc = 0.046 /km. 

The behavior of the phase error for both cases can be observed in Figure 5.4 for two range 
distances using N = 2 entangled NOON states. The phase estimation error depends on the value 
of 4> and it never reaches the Heisenberg limit (the vertical line at S(f> = 0.5). Furthermore, as the 
range increases, the minimum phase error moves away from the Heisenberg limit. 

The behavior of the phase estimation error when one changes the atmospheric attenuation 
by 1 dB/km is shown in Figure 5.5. As can be observed, the effect of this change on the phase 
estimation error is substantial. 

But perhaps more dramatic, it can be observed that the phase estimation error diverges when 
the sine function in the denominator of Equation (5.26) is close to zero. The number of divergent 
points in a range of values of <p depends on the number of photons N. Figure 5.6 compares the case 
of 2, 4, and 8 photons. The reason for these divergences is not physical, but simple a consequence 
that the approximation used for the phase estimation error: 


84 > & 


A Ap 

i d(Ap) , 
1 dtj> 1 


AT | sin A7(0 — 0 O )| 


(5.35) 
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R=100 km R=300 km 

7 = 0.0016 / km 

5<p Ac 



R=5 km 


jj . x = 0 . 046 1 km 

5<p A ' c 



Figure 5.4: Phase estimation error as a function of the phase 0 using N — 2 NOON states. Notice that, 
due to atmospheric attenuation, 8(p never reaches the Heisenberg limit. 

X r = 0 . 0016 / km x c ~ 0.016/fcw 



Figure 5.5: A change of 1 dB/km in the atmospheric attenuation constant produces a significant effect 
on the minimum value achievable by the phase estimation error. 


does not hold well in the regions where the sine function in the denominator becomes zero. That is, 
the approximation does not hold for those points near: 

N (0 — 0o) ~ 0, 7i, In , ..., (!n)n (5.36) 


where n is an integer number. 

In any event, one can see that the phase estimation error is minimal when the sine function 
in the denominator is 1. That is: 


86 


(e) 



(5.37) 
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N=2 N=4 N=8 



Figure 5.6: Phase estimation error for 2, 4, and 8 photons at 10 km range and 300 km visibility. The 
higher the number of photons, the more divergent points are found. 


where we have ignored the contribution from the second arm, assuming it is in a well controlled 
environment inside the sensing apparatus. 


R=100 km R=300 km 

dtp 



X c = 0.0016/ km 

X c = 0.016 /km 



Figure 5.7: Phase estimation error for non-entangled states. 


As shown in Figure 5.7, a similar behavior is observed in the case of quantum interferometry 
using separable states. In this case, however, the minimum value of the phase estimation error is 
given by: 





1 

sfN V 2 “l 


(5.38) 


Notice that in the limit of no attenuation, these last two expressions reduce to the standard 
quantum limit and the Heisenberg limit, respectively. 
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Now let us assume that somehow we reach the minimum phase estimation error. The range of 
the detection system for the entangled (‘FQ and non-entangled (T )le ) cases are shown in Figure 5.8. 
As expected, the entangled system is able to reach the Heisenberg limit (HL) at very small distances, 
while the separable state merely reaches the standard quantum limit (SQL). For an attenuation 
corresponding to a visibility of 300 m the estimation process using entangled states exceeds the 
standard quantum limit up to about 350 km. On the other hand, supersensitivity is restricted to 
about 12 km when the visibility is close to 30 km. 



Figure 5.8: Phase estimation error as a function of the range to the target for a N — 2 entangled Tj, and 
non-entangled sensors. The extent of the supersensitivity regime, where HL < 8cf> < S QL, depends 

on the atmospheric attenuation constant. 

It is important to observe that, for large ranges the sensitivity of the entangled system becomes 
worse than that for the non-entangled system. This is observed to happen at about 550 km in Figure 
5.8 for the scenario corresponding to 300 m visibility, and at about 19 km for the scenario with a 30 
m visibility. That is, for the high visibility case, the use of entanglement brings an improvement in 
phase measurement up to about 550 km. 

In Figure 5.9 we can observe the performance of the system when we use a larger number of 
entangled states in the NOON state. The case for N = 10 is shown. It can be observed that in the 
high visibility scenario, supersensitivity is achieved up to about 180 km, and at about 7 km for the 
low visibility scenario. 

Figure 5.10 compares the phase estimation error for N = 2 and N = 10. As expected, the 
phase estimation error is reduced as N increases. However, the performance of the highly entangled 
states degrades more quickly than for the simpler entangled state. 

Finally, it is worth comparing the performance of quantum interferometry in the radar and 
visible wavelengths. In Figure 5.11 we show the phase estimation error for an attenuation coefficient 
of 0.87 dB/km, which corresponds to X = 550 nm light traversing average atmospheric conditions. 
It can be observed that the phase estimation error is severely compromised by the atmospheric 
attenuation, and the range of super-sensitivity is of about 0.6 km. 
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s 4>,:,, % = 0.0016 / km <%„, Xc = 0 0461 km 



Figure 5.9: Phase estimation error as a function of the range to the target for a N = 10 entangled 
and non-entangled T )J(? sensors. The minimum value of the phase estimation error is smaller than for the 
N = 2 case, but the range of super-sensitivity is also smaller than for the N = 2 case. 



Figure 5.10: Phase estimation error for N = 2 and N = 10. 


All these effects are a consequence of atmospheric attenuation. If the quantum interferometer 
systems could operate in a perfect vacuum for which Xc = 0, then the sensors using NOON states 
would operate at their Heisenberg limit regardless of the range to the target (and the sensors operating 
with non-entangled states would operate at their standard quantum limit at all values of R). 

So far we have been talking about phase estimation, but these results can be easily generalized 
for the case of range-to-target estimation, the natural function of a radar system. The range-to-target 
estimation error can be approximated by: 


SR sa 8<p (5.39) 

lit 

where X is the wavelength of the photons. Therefore, all our observations for the behavior of 8<fi hold 
the same for SR. 
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Figure 5.11: Atmospheric attenuation for the phase estimation error using visible light (k = 550 nm) 
instead of microwaves. 


One may find it surprising that if we use k = 3.2 cm, then the range error in the supersensitivity 
regime will also be in the order of cm. It is not expected that a quantum radar will operate at this 
sensitivity when detecting incoming targets. Indeed, this analysis has considered the target as a 
perfectly reflective and static object, which allows for a very sensitive measurement. In a sense, the 
scenario discussed so far resembles the situation with the Lunar Laser Ranging Experiment, which 
is able to measure the distance from the Earth to the Moon with milimetric accuracy [109]. 

The fact that the super-sensitivity regime may reach up to 350 km suggests that long-range 
phase estimation maybe feasible using quantum interferometry with NOON states. However, many 
issues remain problematic. For instance, we have assumed that the detectors have perfect sensitivity. 
In reality this is not the case, and the imperfect performance of the sensors involved will diminish 
the supersensitivity regime achievable by the system. In addition, the expressions derived present 
periodic divergences. 

As a consequence, these observations make evident that long-range phase estimation using 
basic quantum interferometry with NOON states may be severely restricted due to the attenuating 
effects of the atmosphere. However, this only means that NOON states alone are not enough to 
build a robust interferometric quantum radar for practical applications. 

5.2.5 ADAPTIVE OPTICS CORRECTION 

As we have seen, because of atmospheric attenuation, the use of NOON states is not enough to 
guarantee supersensitivity on a quantum radar. To circumvent this problem, James Smith from 
the US Naval Research Laboratory developed an adaptive optics correction method that allows 
supersensitivity over ranges up to 5000 km and for atmospheres with significant variation in their 
electromagnetic properties [97]. 
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The proposed atmospheric correction method is inspired by the guide stars correction tech¬ 
nique, which is widely used in adaptive optics to improve the performance of telescopes [101]. In 
general, a turbulent atmosphere is not a homogeneous medium. This causes fluctuations on electro¬ 
magnetic properties such as the optical refractive index. These variations are responsible for wavefront 
distortions that produce loss mechanisms such as beam spreading, beam wandering, and intensity 
fluctuations (scintillations). In astronomy, the correction of images of dim stars is accomplished using 
a nearby bright star as a reference to determine the shape of the distorted optical wavefronts. It is 
assumed that the light from both sources pass through the same portion of the atmosphere. The 
telescope uses a wavefront sensor to measure the distortions that the atmosphere has introduced to 
the incoming stellar light. Subsequently, a computer determines the optimal shape of an mirror that 
can correct for these distortions, and a flexible mirror is shaped accordingly. 

This atmospheric correction technique can also be used to improve the performance of RF 
and laser radar [110, 111]. In the former case, the “guide stars” could be implemented using small 
electronic devices that produce electromagnetic radiation of fixed frequency and amplitude. These 
devices can be fired by artillery and could provide a mean to measure the refraction, attenuation, and 
other electromagnetic properties of the atmosphere along the line-of-sight path between the source 
and the detector. These values can be extrapolated over time and space to provide an approximate 
electromagnetic map of the atmosphere that surrounds the radar and the target. In addition, the 
adaptive optics system may be fed with historical information as well as with data measured from 
electromagnetic sources of opportunity. 

The adaptive optics technique in a quantum radar introduces a phase shift 0„ o that can be 
controlled by the user. Then, neglecting the contribution from the small interferometer arm L i, the 
corrected phase estimation error is given by: 


<50 


(c) 


~ 2 + sin2 N((p + * ao ~ W 


N | sin IV(0 + (p ao — 0q C) ) I 
where the corrected dispersion phase shift is: 

0n } ^ - ] 1 {C) l 2 


the corrected transmittance of the medium is: 


» 


„-x (c) l 2 


(5.40) 


(5.41) 


(5.42) 


and the parameters r] (<y> and / are the corrected refractive and attenuation constants of the medium, 

respectively. 

As made evident in the analysis of atmospheric quantum interferometry, the phase estimation 
error diverges when the sine function in the denominator becomes zero. Therefore, the adaptive 
optics technique also requires for the user to implement a controlled phase delay 0 ao in a way such 
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that <50 (< 4 is minimum when we estimate the range to the target Li. That is, cj) ao has to be chosen 
in a way such that the sine function in the denominator becomes 1 or — 1: 




t r 3 tt (2n + \)n 


Therefore: 




1 1 

+ - 


N y 2(a\ c) ) N 2 


and similarly for the non-entangled case: 
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V~N Y 2 a[ c) 2 


(5.43) 

(5.44) 

(5.45) 


Unfortunately, the adaptive optics technique for quantum interferometry requires of estimated 
values for several parameters, including the range to target Z. 2 - In a sense, the adaptive optics 
technique for quantum interferometry implements an information fusion and filtering process that 
improves the values of the estimated variables. This means that adaptive-optics interferometric 
quantum radar may be better suited to discriminate or identify targets that have already been detected 
by other means. 

In any event, it is expected that the use of adaptive optics in an interferometric quantum radar 
will be able to increase the supersensitivity regime to well over 1000 km [97]. 


5.3 QUANTUM ILLUMINATION 

Quantum Illumination is a revolutionary photonic standoff quantum sensing technique developed 
by Seth Lloyd at MIT [67]. This technique enhances the sensitivity of photodetection of light 
in noisy and lossy environments. By design, quantum illumination is not restricted to any specific 
frequency, so in principle, it could work in the radar X-band regime. 

The setup is similar to the case of single-photon and entangled-photons quantum radar 
discussed above. We assume a target that is illuminated by a quantum illumination device. The 
goal is to detect the presence of the target even in a noisy and lossy environment. As we will see, 
entanglement increases the sensitivity of the detection system. Quite amazingly, the enhancement 
provided by quantum illumination using entangled photons is only observed in a noisy and lossy 
environment. 

Quantum illumination has a different detection strategy than an interferometric quantum 
radar. In the present case, one does not perform phase measurements, but simple photon detection 
counts are enough. Nevertheless, this is not an easy task. As we will discuss later, although the 
measurement of entangled photons is physically realizable, it is a very challenging procedure. 

The interaction between the signal photons and the target can be modeled as a beam splitter 
with small reflectivity r]. Noise is injected into the system with average photon number b into each 
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optical mode (frequency or polarization). Although it is not necessary, we can assume that b <5C 1 
to simplify the mathematical formalism, but the general results remain valid for arbitrary b. The 
b 1 case corresponds to the scenario where the thermal radiation is substantially below the signal 
photon energies. 

If the photodetector is characterized by a bandwidth W and a temporal detection window T , 
then, in principle, the device can distinguish between d different radiation modes: 

d^WxT (5.46) 

Furthermore, we assume that T is small enough so that at most one noise photon is observed per 
detection event, that is: 

d x b « 1 (5.47) 

We consider two modes of operation, when the signal photons are not entangled and when 
they are entangled. These two modes will be analyzed and compared in what follows. 

5.3.1 NON-ENTANGLED PHOTONS 

Let us first assume the case of non-entangled signal photons. Two cases are possible: (1) the target 
is not within range, so it will be undetectable by the system (but it is possible that the system will 
record a false positive), and (2) the target is within range, and there is some probability that the 
system will succeed in detecting it. To analyze the performance of quantum illumination we need to 
compute the probabilities for all these events. 

Case (1): No Target Within Range 

If the target is not within range then the signal photon will be lost. The only photons detected by 
the system are noise photons which can be described by: 

Pi = Pb ® Pb ® ■■■ ® Pb (5-48) 

v -v-' 

d states 

Indeed, there is one possible noise photon state with b photons on average {pb) for each of the d 
possible optical modes. We can rewrite pi as: 


Pi = ® Pb ’ 


k= 1 


where: 

Pb ] = MW 


(5-49) 


(5.50) 


is the density operator for a single noise photon in the mode k. 
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Because the average number of photons per detection event is small, db 1, we can approx¬ 
imate p i as: 


Pi *(l-db)\0)(0\ + bJ2\k)(k\ 


(5.51) 


k= l 


where | k) is a single-photon state in the mode k and the vacuum state is represented by |0). That is, 
there is a probability (1 — db) that we will not measure any signal, and a probability b that we will 
measure one out of d possible noise photons. 

The accuracy of the approximation can be computed through the norm distance, which results 
in: 


approx ,0, II ~ i 

and therefore, the approximate expression is correct as long as we remain in the lowest order of b. 


p\ rr "" - pb® Pb® ■■■ <8> Pb\\ ~ d 2 b- + OQf) (5.52) 


Case (2): Target Within Range 

If the target is within range of the quantum illumination system, then there is a probability i] that 
we will measure a thermalized signal photon p, and a probability (1 — p) that we will just measure 
noise pi.That is: 

Pi = (1 - p)p\ + pp (5.53) 

where p\ is given by Equation (5.51). However, it can be shown that: 

*7 IIP ~ Pll i ^ pdb + 0{rjb 2 ) (5.54) 


where p is the signal photon state. In the regimen of interest we can ignore terms O(iib) and as a 
consequence: 


P ^ P 


(5.55) 


which leads to: 


Pi ^ (1 - P)P 1 + PP 


(5.56) 


Detection Probability 

We have a probabilistic problem with two possible hypothesis: (TR) there is a target within range 
or (NTR) there is no target within range. To discriminate between these two hypothesis we need to 
repeat the probabilistic experiment several times. That is, we need to repeatedly send signal photons 
towards a region of space to determine if there is a target. 

Using the expressions for p\ and p 2 we can determine the probabilities of detecting a target 
within range, measuring a false positive (detecting a target that is not there), or the probability that 
no target is present in the region of interest. 

Thus, for instance, the probability of a false positive is the probability of positive detection “+” 
of a target that is not within range, NT R. This is obtained from p\ by computing the probability of 
measuring one of the d noise photons in the mode k. This probability is given by: 


p(+\NTR)=b 


(5.57) 
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On the other hand, the probability of getting a negative detection ” (that is, the probability of 
successfully concluding that there is no target within range) is given by: 

p{-\NTR) = \-b (5.58) 


A similar analysis follows for the case where there is a target within range, T R, and the entire 


set of probabilities are: 

p(-\NTR) = 
p(+\NT R) = 
p(-\TR) = 
p(+\T R) = 


1 - b (5.59) 

b 

(1-feXl-ij) 
b( 1 - v) + 9 


which clearly satisfy the probability normalization conditions: 

NTR : p(-\NTR) + p(+\NTR) = 1 (5.60) 

TR : p(-\T R) + p(+\T K) = \ 


The number of times we need to repeat the probabilistic experiment is the number of photons 
that we need to send towards a region of space to determine if a target is within range. This number 
depends on the signal-to-noise ratio (SNR) of the quantum illumination system and can be estimated 
to be: 


SNR 


P(+\TR ) 

P(+\NTR) 
b( 1 — rj) + 77 

b 

n 

1 - r) + - 

O ( rj/b ) 


(5.61) 


Depending on the value of SNR, two regimes of interest have to be considered: (H) a High 
SNR regime in which SNR> 1, and (L) a low SNR regime in which SNR< 1. For the first case, 
the number of photons required to discriminate between the two possible hypotheses is given by 
Nh and is approximated by: 

N h ~ 0(l/r,) (5.62) 

On the other hand, for the low SNR regime, the number of photons necessary to discriminate 
between the hypotheses is denoted by AT and given by: 

AT ~ OiSb/p 2 ) (5.63) 


As expected, Nh is smaller than AT: 

AT « Otfb/p 2 ) « 0(b/p) x 0(1/??) « Q(b/p) x N H > N H (5.64) 


in the low SNR regime (b > ??). 





5.3.2 ENTANGLED PHOTONS 
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Let us now consider quantum illumination using entangled photons. As usual, we need to generate 
entangled signal and ancilla photon states. The signal photon (S) is sent towards a target and the 
ancilla photon (A) is kept within the radar system. Both photons are entangled in their momentum 
degrees of freedom over all possible modes d: 

1 c> 

\^ e) = ^Y J \Vs\k) A (5.65) 

k =l 

Recall that \k) represents single photon states in the mode k. As a consequence, represents a 

2 photon system. The associated density matrix for the entangled state is: 

p e = np)(«)M(ip| (5.66) 

As in the case for non-entangled photons, we have two cases to consider: (1) when a target is 
not within range and (2) when the target is within range. 


Case (1): No Target Within Range 

If there is no target within range, then the system will just be able to detect noise photons. The 
signal state is considered lost and then it can be approximated by p\ as described in the case for 
non-entangled photons. At the same time, the ancilla state goes to the completely mixed state: 


P\ 


+ 0(b 2 ) 

d 


(5.67) 


(l-db)\0) s (0\+bJ2\Vs(k\ 


k= 1 


I ^ + 0(b 2 ) 
d 


Ia 


[(1 - Jh)|0) 5 (0| + bis] ® ^ + Oib 2 ) 

d 

and Is and I a are the identity operators on the single photon Hilbert space for the system and ancilla 
photons, respectively: 


k =1 


(5.68) 


for X = S, A. 


Case (2): Target Within Range 

If the target is within range, the system will detect the original entangled state with probability 
rj, and it will detect a noise photon p\ with probability (1 — rj). As before, the signal state will be 
thermalized, but it can be represented by the original state if we are kept within terms of the order 
of 0(b, b 2 ). That is: 


/o| = (1 — rj)p[ + r\p e + 0(r]b, rjb 2 ) 


(5.69) 
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where once again we have replaced the thermalized density operator p e by p e , an operation that can 
be done in the approximation regime under consideration. 

Detection Probability 

As in the case with non-entangled photons, we need to discriminate between two hypothesis. The 
associated probabilities can be computed from p e x and p|. Then, the probability of positive detection 
“+” of a target that is not within range, NT R, is obtained from Pj and is given by: 


p e (+\NTR) = - 
d 


(5.70) 


Similarly, the probability for a negative detection ” (that is, correctly concluding that there is no 
target within range) is: 



These two equations make clear the effect of entanglement, as they are the same for the case of 
non-entangled photons but making the replacement: 



(5.72) 


In other words, the noise b is effectively reduced by a factor d. 


A similar conclusion follows for the case where the target is within range. The entire set of 
probabilities is then given by: 


p e (-\NTR) 



(5.73) 


p e (+\NTR) 


b 

d 



which clearly satisfy the probability normalization conditions: 


NTR : p e (—\NTR) + p e (+\NTR) = 1 
TR: p e (-\TR) + p e (+\TR) = l 


(5.74) 


For the entangled-photons case, the signal-to-noise ratio is given by: 


SNR e Q(r]/b x d) 


(5.75) 


which reflects the fact that the entanglement has increased the SNR by a factor of d. 
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Once again, two regimes appear, a (H) high SNR regime where SNR> 1 and a (L) low SNR 
regime in which SNR< 1. The number of photons necessary to determine the presence of a target 
in the high SNR regime is given by: 

N e H « 0(l/r,,) (5.76) 

On the other hand, for the low SNR regime, the number of photons necessary to discriminate if a 
target is within range is: 

N e L « 0(Sb/drj 2 ) (5.77) 

5.3.3 SENSITIVITY COMPARISON 

All the relevant results for quantum illumination are condensed in Table 5.2. We can conclude that 
the use of entangled states increases the signal-to-noise ratio by a factor d: 

& d > 1 =>■ SNR e > SNR (5.78) 

SNR 

The advantage of a higher SNR is its effect on the number of photons necessary to elucidate 
the presence or absence of a target in a given region of space. As seen in the table, for the high SNR 
regime, entanglement does not improve the performance of the system. It is only in the low SNR 
regime where entanglement provides an advantage over single-photon pulses. Indeed: 

AT 1 

-rr « 7 < 1 =» N e L < N l (5.79) 

N l d 


Table 5.2: Sensitivity of Quantum Illumination for Entangled 
and Non- Entangled States. 


Non-Entangled Photons 

Entangled Photons 1 

p(-\NTR) 

1 -b 

1 — b/d 

P(+\NTR) 

b 

b/d 

P(-\TR ) 

(1 — b)(\ — rf) 

(l-b/d)(l-ri) 

P(+\TR) 

b( 1 — rj) + rj 

(1 — r\)b/d + 11 

SNR 

0(r)/b ) 

0(r\d/b ) 

High SNR 

il/b > 1 

rjd/b > 1 

N h 

0(1/1 7) 

0( 1/ri) 

Low SNR 

r\/b <\ 

rjd/b < 1 

Nl 

0(8b/i 7 2 ) 

0(8b/rj 2 d) 


Let us recall that an ebit is one unit of bipartite entanglement. In other words, an ebit is the 
amount of entanglement that characterizes a maximally entangled state made of 2 qubits. Thus, 
a 2-qubit Bell state contains 1 ebit of entanglement. And in general, a maximally entangled state 
made of two systems of dimension d has m = log 2 d ebits of bipartite entanglement. Therefore, 
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quantum illumination in the low SNR regime increases the sensitivity of the detector by a factor 2 m . 
That is, considering entanglement as a resource, quantum illumination offers an exponentially large 
improvement on the sensitivity of the detection system. 

It is also important to observe that, within the approximation considered, the reflected signal 
photon can be considered as fully thermalized. In other words, any possible level of entanglement 
between the signal and ancilla photons is annihilated by thermal effects. But even so, quantum 
illumination is able to extract the correlation between these states to enhance detection. By any 
means, this is a rather surprising and unexpected result. 

As previously mentioned, the treatment presented here is limited to the low-temperature 
regime which corresponds to a small number of noise photons per mode (b 1). However, quantum 
illumination has been shown to work even in the high-temperature case [67]. The final result is the 
same: sensitive performance is enhanced in the low SNR limit by a term exponentially large in the 
number of ebits of entanglement used for the photon states. 

Finally, it is worth noticing that the calculation of detection probabilities presented here 
was limited to the use of optimal single-shot measurements [44] . That is, we just computed the 
projection of the desired outcome over the quantum state. However, a formal treatment has to 
involve the Chernoff Quantum Bound which determines the optimal strategy that minimizes the 
total probability of error during the discrimination of two different quantum states when one can 
access an asymptotically large number of copies [4, 12]. A formal estimation of the probabilities 
using the Chernoff bound can be found in Lloyd’s paper and leads to the same results presented 
here [67]. 


5.3.4 GAUSSIAN STATES 


The study of quantum illumination can be extended for the case of entangled Gaussian or coherent 
states [48, 94]. In the number-ket representation, a two-mode (system S + ancilla A) monochromatic 
Gaussian states takes the form: 


\'f)sA = 


Ei 

n =0 




(N s + 1)' I+1 


\n)s\n)A 


(5.80) 


where Ns is the average number of photons per mode [73, 92]. The above expression is given in 
terms of the number-ket representation. That is, the states \j) and \i) represent states of j and i 
photons, respectively, with some common momentum k. 

It can be shown that by making an optimal joint measurement on the entangled Gaussian 
state, quantum illumination offers a 6 dB improvement in the error probability over a coherent state 
system in the lossy, noisy, low-brightness scenario [48, 94]. 

In addition, it has been found that quantum illumination with entangled Gaussian states 
brings increased performance in the angular resolution for the reliable discrimination between one 
or two possible targets sustaining a small angle as viewed by the transmitter [37]. 
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5.3.5 ENTANGLED MEASUREMENTS 

Arguably, the most challenging aspect of implementing a quantum illumination radar system is the 
ability to perform detection of entangled states. These types of measurements have been successfully 
performed in the laboratory, but they are difficult to accomplish [18, 28]. Indeed, the detector has to 
verify that the signal and ancilla photons are correlated in their frequency and time of arrival [67]. 

The correlation on their time of arrival may be difficult if the range to the target is unknown. 

In such a case, the quantum illumination system has to delay the ancilla photon across a large 
number of times until it can find a match. This added complexity suggests that perhaps quantum 
illumination is best suited to discriminate the number of possible targets once a positive detection 
has been established using traditional radar or other sensing systems. 

Furthermore, the detector has to verify the energy correlation without distinguishing between 
the different possible modes k: 

col + = co (5.81) 

where co is the frequency of the original photon used to create the entangled pair through parametric 
down - co nversio n. 

In spite of these challenges, a feasible quantum illumination detector has been recently pro¬ 
posed in the scientific literature [38, 39]. 

5.4 QUANTUM RADAR JAMMING 

As we discussed in Chapter 4, active radar countermeasures include the use of jamming to hide the 
presence of a target. In this regard, it is important to mention that quantum radar appears to be 
more robust against the effects of a jamming signal. Indeed, let us recall that both quantum radar 
proposals, the interferometric quantum radar and quantum illumination, take advantage of the same 
principle to achieve higher performance. That is, these systems perform collective measurements 
after highly correlated states have been injected into the system. These correlations can be exploited 
to prevent signal photons to get lost in noise or jamming signals. 

Let us consider, for example, the case of quantum illumination using non-entangled photons. 
The signal-to-jam ratio will be approximately given by: 

SJR ne « 0(rj/j) (5.82) 

where i] is the reflectivity of the target and j is the average number of photons injected by the jammer 
into the system. 

On the other hand, if we use entangled photons, the signal-to-jam ratio is now given by: 

SJR e « 0(T]/j X 2'") (5.83) 

where m is the number of ebits of entanglement in the photon states. 

As a consequence, entangled quantum illumination is less likely to be jammed than a quantum 
radar operating with non-entangled photons. That is, in terms of entanglement as a resource, a 
quantum radar is exponentially less likely to be jammed than a classical radar. 
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5.5 PHYSICAL REALIZATION OF A QUANTUM RADAR 

To date, quantum radar remains an intriguing theoretical proposal. The hardware that enables quan¬ 
tum radar transmission and reception remains to be implemented. While the generation and detec¬ 
tion of entangled photons in the visible regime is done quite easily with today’s technology, it is not 
yet possible to perform these operations with photons in the microwave range. 

However, some ideas have been proposed towards the generation and detection of photons in 
the radar region. Some of these will be briefly discussed in what follows. 


5.5.1 TRANSMITTER 


The key element in the quantum radar transmitter is the generation of entangled photons in the 
microwave regime. Parametric down-conversion is the most widely known technique to entangle 
photons in the visible regime [84]. This entanglement generation technique uses a non-linear crystal 
that splits an incoming photon into two entangled photons. 

A less known entanglement generation strategy is the use of semiconductor nanostructures. 
In particular, the decay of bi-excitonic states in quantum dots involve inter-band, transitions be¬ 
tween valence and conduction bands that generate entangled photons with frequencies in the visible 
range [5, 99]. 

However, a similar technique has been proposed for the generation of entangled photons 
in the microwave regime [27]. This strategy produces entangled photons produced by intra-band 
transitions of conduction band electrons. Photon-assisted tunneling experiments have shown that 
these transitions couple to microwave photons [108]. As a consequence, microwave photons are 
produced from spontaneous downward transitions between single-particle levels in a quantum dot. 

This entanglement generation scheme can be accomplished by using an array of four quantum 
dots connected to two electron reservoirs and inserted inside a cylindrical microwave resonator as 
shown in Figure 5.12 [27]. For convenience, the dots are labeled (T) top, (B) bottom, (L) left, and 
(R) right. The two lateral dots, L and R, are used to provide unique initial and final states for an 
electron in the conduction band. On the other hand, the T and B dots provide two decay paths. 

The electron transitions are shown in Figure 5.13. The electron begins in the quantum dot 
L and tunnels into excited states of the dots T and B. The state at this stage can be assumed to be 
symmetric and given by: 


\xfr) = — (|7i) + | Bi)) 
V2 


(5.84) 


That is, the electron can be in the excited state ofT or B with equal probability. Then, the electron 
state decays into the ground states of T or B and in the process emits two photons. If the electron 
decays from T, then both photons have left-handed circular polarization states. On the other hand, if 
the electron decays from B, then both photons have right-handed circular polarization states. Then, 
the quantum state at this stage can be described as: 
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(5.85) 
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Figure 5.12: Arrangement of four quantum dots for a microwave photon entangler [27]. 
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Figure 5.13: Electron transitions in the quantum dots that generate a pair of entangled microwave 
photons [27]. 

which can be rewritten as: 

where the electron states are: 

and the photonic states are: 

The coupling with the dot R serves to disentangle the photons from the electron states. Furthermore, 
the antisymmetric state \4 r ~) is suppressed due to destructive interference. Then, the state of the 


ivp> = 4=(i^+)i<i> + }+i^_)i<i>_)) 

V2 

(5.86) 

W±) = 4=(|7b)±|Bo» 

V2 

(5.87) 

|4>±) = 4=(| + +)±|--)) 

(5.88) 
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photons is found by finding the projection of the state | 'k) with | i/r+). That is, the outgoing photons 
are described by: 


W+l*> ex |«D+) = 


1 


(! + +) +1 


-» 


(5.89) 


which is a state with maximal entanglement. 


5.5.2 RECEIVER 



Figure 5.14: Energy levels and decay paths of proposed microwave sensing metamaterial [88, 89]. The 
device absorbs a photon of frequency co and has a decay rate T to the long-lived stable ground state |g). 

Photon detection is a well known technology widely used in the visible and near-visible 
regimes [19, 25, 83, 87]. However, the design and development of single photon detectors in the 
microwave regime present a large number of theoretical and experimental challenges [88, 89]. Some 
of these difficulties include the unavailability of cryogenic linear amplifiers capable of working with 
low photon number counts, as well as the problem of synchronizing the detection process with the 
arrival of the measured field. 

However, a novel sensing technique has recently been proposed that circumvents most of the 
problems associated to the task of detecting single photons in the microwave regime [88, 89]. The 
idea rests in the design of a novel metamaterial 4 . The proposed metamaterial is designed to absorb 
microwave photons and can be used to perform single-shot photodetection. 

The proposed design of the system consists of a linear waveguide coupled to a set of three-level 
photo-absorbing particles in arbitrary positions. The three energy levels of the absorber are depicted 
in Figure 5.14. It can be shown that these absorbers are able to capture a microwave photon from the 
initial state |0) and transition to the excited state 11). Subsequently, the state decays into a long-lived 
stable ground state |g) with a decay rate T. The photon counting then boils down to a counting of 

^Metamaterials are artificially engineered materials designed to have specific physical properties. Furthermore, the mechanical and 
electromagnetic properties of a metamaterial are a consequence of its molecular structure rather than its composition. 
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those absorbers found in the ground state |g). Therefore, this detector resembles photographic film 
in the sense that once a photon has been absorbed by the metamaterial, the device is changed into 
a stable and mesoscopically distinguishable state. 

One potential problem with the proposed design is its small operational bandwidth, which 
has to be small in comparison to the time it takes to absorb a photon. That is: 

A co < — (5.90) 

However, this may not be an important issue in its application towards a quantum radar receiver. 

Another potential problem is that if the decay process 1 —> 0 takes place, then there will be 
incoming photons that will not be counted. However, recent experiments appear to suggest that 
this process happens with a rate of a few MHz, and therefore it will affect long wavepackets [45]. 
It remains to be seen if this issue could present a serious limitation to the use of these devices for 
quantum radar applications in the X-band. 

5.6 SUMMARY 

Quantum radar is a standoff sensor that could be used to beat the standard quantum limit of current 
sensing devices. Two major designs have been proposed: interferometric quantum radar and quantum 
illumination. The hardware to enable quantum radar is still in the design process, but some work has 
been done towards the development of entanglement generators and photodetectors able to operate 
in the radar region. 

Unfortunately, to date it is very difficult to predict the performance of a quantum radar in an 
operational environment. First, most quantum radars proposals presume a perfectly reflecting target 
oriented in the specular direction. Second, we cannot use the standard definition of the radar cross 
section. 

Indeed, regardless of the underlying system or states used, the important point of quantum 
radars is that, as their name suggests, they operate on the quantum regime. That is, they use a relatively 
small number of photons. As a consequence, theoretical consistency requires that the electromagnetic 
fields have to be described as quantum fields. Thus, the radar pulse is not described by classical 
electrodynamics, but by quantum electrodynamics. Therefore, just as the radar cross section was 
defined through Maxwell equations for classical fields, it is necessary to define a quantum radar cross 
section ctq based on quantum electrodynamics. 
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CHAPTER 6 

Quantum Radar Cross Section 


So far our discussion on quantum radar has ignored the geometry and composition of the target. In 
practice, of course, targets will present intricate geometries that will reflect the incoming photons 
on a complex pattern. Within the realm of classical radar theory, the radar cross section ac is used 
to determine the “radar visibility” of a specific target. 

We recall from Chapter 4 that ac is defined as the ratio between the power reflected back 
to the receiver per unit solid angle and the incident power density. As we discussed, ac emerges as 
a critical concept in classical radar theory, as it offers an objective measure of performance for the 
radar system and the stealth capabilities of modern weapon platforms. 

As quantum radars emit a handful of photons at the time, the radar-target interaction in this 
regime is described through photon-atom scattering processes governed by the laws of quantum 
electrodynamics. As such, it is theoretically inconsistent to use the same ac to characterize the 
visibility of a target illuminated with a quantum radar. As a consequence, there is a need to develop 
the concept of a quantum radar cross section aQ to objectively measure the “quantum radar visibility” 
of a specified target. That is, we need to define erg to analyze the performance of quantum radars 
in the scenario where the targets are not perfectly reflective objects and the radar signal is made of 
a handful of photons [61]. 

6.1 DESIRED FEATURES OF a Q 

There are several conceptual challenges that make it difficult to define erg. For instance, energy in 
the classical regime is proportional to the amplitude. But in the quantum realm the energy of a 
single photon is proportional to its frequency. Furthermore, as we discussed in Chapter 2, photons 
are nonlocalizable. 

However, we can attempt to establish the basic properties that are desired in a conceptually 
robust definition of erg. In simple analogy to the classical radar cross section, the quantum radar 
cross section should have these properties: 

Operational Meaning: The reason to define a g in the first place is to have an objective measure 

of the quantum radar visibility of a target. 


Energy Conservation: As with ac, aQ should entail energy conservation in the optical regime 
when absorption effects are ignored. 
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Strong Dependencies: Similarly, to cre > it is desirable that oq strongly depends on properties of the 
target: geometry (absolute and relative size, shape and orientation), as well as its composition 
(material properties). 

Weak Dependencies: In the same manner, it is desired that oq is approximately independent on 
the properties of the radar system. That is, oq should depend very weakly on the strength, 
architecture, physical implementation, and range of the radar system. 

Multiplatform Comparison: To better understand the advantages and disadvantages of quantum 
radars, it would be desirable to be able to directly compare oq and oq. 

Asymptotic Behavior: In the large photon limit, the quantum realm gives way to classical physics 
and oq should be proportional to oq: 

lim Oq OC Oq (6.1) 

n y —>OQ 


6.2 INCIDENT AND SCATTERED QUANTUM FIELDS 

As we discussed in Chapter 3, it is possible to describe all the physical properties of specular reflections 
exclusively using QED processes of absorption and emission of photons by the atoms in a mirror [62]. 
Following an interferometric analysis, we found that the intensity measured by a detector after a 
photon is reflected by N atoms is given by: 


(l s its, t)) 
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( 6 . 2 ) 


where the photon wave functions are: 

'U, /) (A Rj, t) = - A Rj/c) e -0'®+ r / 2 )6-A«,/c) 

Y Ar id 
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and T represents the inverse of the lifetime of the excited state of the atom: 


(6.3) 


(6.4) 
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(6.5) 


In the above expressions &> is the frequency of the incoming photon, and i] is the angle between the 
electric dipole moment of the atom J1 and Ar ^ = r; — is the distance between the i th atom and 
the detector. In the context of quantum radar, A Rj is the total interferometric distance from the 
radar transmitter to the target, and from the target to the radar receiver. 








6.2. INCIDENT AND SCATTERED QUANTUM FIELDS 95 
For monostatic radars the source and the detector are at the same spatial location: 

r s = r d ( 6 - 6 ) 

and therefore 


A Ri = A r si + Arid 

= |r,s -r;| + |r; -iy| 
= 2|r,--iy| 


(6.7) 


where r s , rj, and r,- are the positions of the source, the detector, and the i th atom, respectively. 

From an operational point of view, the definition of oq has to be able to provide an objective 
measure of the quantum radar visibility of the target. Therefore, oq has to depend explicitly on 
( I s {r, 6, (f>, t)), the scattered intensity at the receiver: 


cr q oc {I s (r, 0, <p, t)) 


because in the large photon limit: 


lim ( I s (r,9,4>,t )) oc |E J@r | 
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(6.9) 


On the other hand, the resulting average intensity incident over all the N atoms that make 
the target is (assuming the target behaves as a photodetector): 
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where we have presumed that in the region of interest: 

A r si = |r. s -Ti\& R Vi (6.11) 

It is important to note that we can make this approximation because the exponential terms in the 
equation of (/, ( r s , t )) do not form an interference pattern. Indeed, the absolute value is inside the 
summation symbol. In other words, there is only one possible path from the emitter to the target. 

The © step function merely indicates that the signals will not propagate faster than the speed 
of light. With this understanding, we can remove it from the equation to simplify the notation. 
Therefore: 

(h(r s J)) ~ e- r{, - R ' c) ( 6 . 12 ) 
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and as a consequence, the integration over the orthogonal projected area of the target gives: 


L 


(ii(r s ,t)) ds * A ± (e,4>) (ii(t s ,t)) 


(6.13) 


T±(e,<p) 


If we ignore absorption effects, energy conservation in the optical regime requires that all 
the incident energy has to be reflected in some direction. Thus, we can derive expressions similar 
to those obtained in Section 4.9. That is, integrating the position of the detector over a spherical 
surface Sr cI of radius R d : 
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which can be approximated by: 
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6.3 OPERATIONAL DEFINITION OF a Q 

In analogy to gc, it is reasonable to define gq as: 

o B s lim 4nR 2 (6.16) 

(//(r„0) 

and if we assume energy conservation in the optical regime, then oq for a monostatic quantum radar 
can be approximated by: 


gq rs 4jtA±(0, 4>) lim 
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(6.17) 


where the expectation value of the scattered intensity is taken at the receiver [61]. 

Unfortunately, this “simplified” expression for gq remains problematic for the analytic study 
of quantum radar. In general, numerical methods have to be employed to elucidate the behavior of 
the quantum radar cross section. 

It is possible, however, to get an analytical expression for the quantum radar cross section 
of a spherical target. Indeed, for this target is safe to assume that Gg l>he ' ‘ 1 will showcase spherical 
symmetry and it will be independent of the angular variables 0 and <p. Then: 

jf" JJ a Q Phere) sin °Md4> = AitGf here) (6.18) 

In addition, the expression for Gq 1 ’ 1 ’" ‘ * in the optical regime leads to: 
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( 6 . 19 ) 





and as a consequence: 
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(. sphere ) 

a Q 


A± 


( 6 . 20 ) 


which means that, in the case of a spherical target, the classical radar cross section and the quantum 
radar cross section are the same: 


(, sphere ) ( sphere ) 

°Q ~ a C 


( 6 . 21 ) 


In general, for non-spheric targets, the quantum and the classical radar cross sections maybe different. 


6.4 THE QUANTUM RADAR EQUATION 

We are now in position to derive a quantum radar cross section based on the definition of quantum 
radar cross section. From the previous definition for the quantum cross section we have that: 

an = lim 4nR 2 ^-P- (6.22) 

R — >oo (/,-) 

and using the previously obtained expressions for the intensities we get in the large R limit: 

(/,) * (6.23) 
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This last expression resembles the classical radar equation. Indeed, if we make the replacements: 

T f G = 4 t r El (6.24) 

P, Q = (I s )A r 

where P , G is the transmitted power of the quantum radar and Pp is the received power at the 

quantum radar. Then: 

pQ — Pi _ ^ r (6 25) 

r ~ ( 477 ) 2^4 

which for all means and purposes can be considered as the quantum radar equation. 


6.5 SIMULATION OF FOR RECTANGULAR TARGETS 

We have developed computer programs to simulate oq and oq for a finite, flat rectangular plate. 
The geometry of the sensor-target system is the same as the one shown in Chapter 4 and shown 
again in Figure 6.1. As before, we will ignore diffraction and absorption effects. 
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Figure 6.1: Sensor-target geometry. 


Even though it is not computationally efficient, we decided to perform the coding using Math¬ 
ematical M because it allows the rapid prototyping of simulation programs [75]. As a consequence, 
the target size and atomic composition had to be severely limited by the available computational 
resources. 

For instance, a oq simulation of a target made of 10,000 atoms and a single quantum radar 
signal photon took over 4 hours on a 1.6 GHz Intel desktop. Furthermore, multiple incoming radar 
photons contribute an exponentially large number of terms in the simulation. Thus, our computer 
simulations were limited by 5 radar photons and 100,000 atoms. By any means this is a very simple 
radar-target system. Nevertheless, this model offers a useful way to obtain valuable insight on the 
physical properties of oq. 

Finally, it is interesting to mention that the simulation of oq appears to be an ideal task to be 
solved by quantum computers [63] . 
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6.5.1 ANALYTICAL EXPRESSION 

Using the equations for the expectation value of the intensity we obtain: 
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where we have assumed that outside the interference terms, we can approximate: 

Ar ;d = \t d — r,-1 ^ R Vi 

in the amplitudes. 

We can further assume that T & 0, which leads to: 


(6.27) 
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6.5.2 SIDELOBE STRUCTURE 

Figure 6.2 shows the plot of erg vs. 0 for 0 = 0. Similarly, to ere, the highest value of erg is reached 
when the target is oriented in the specular direction at 9 = 0. In addition, a sidelobe structure is 
observed. However, while the sidelobe structure for oq was due to end-region returns, the sidelobe 
structure for erg is a purely quantum mechanical effect due to quantum interference. 
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Plate Size: 



Figure 6.2: Simulation of erg vs. #(</> = 0) for a 2D flat plate. As with oc, the highest value of erg is 
reached when the target is oriented in the specular direction. The sidelobe structure, however, is a purely 
quantum mechanical effect. Diffraction and absorption effects have been ignored. 


6.5.3 oq VS. oq 



Figure 6.3: Comparison of erg vs. cre for a rectangular target. 
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It is possible to compare oq and oq because both have the same units of area and they have 
the same operational definition of ratio of intensities. The result is shown in Figure 6.3. It can be 
observed that oq is approximately equal to oq for most values of the angular variable 0. 

We need to remark that this result is valid only for the rectangular target. In addition, let us 
recall that for a spherical target, oq = oq. However, it remains to be shown what is the behavior of 
oq in relation to oq for singly curved surfaces. This is something we will explore in the future. 

We also observe that the values of oq and oq are different towards the region of large scattering 
angles. This, however, may be due to the error in the approximation used in the equation for oq. 
Indeed, let us recall that the physical optics approximation tends to fail at large scattering angles. 
For example, a similar behavior is observed in the classical approximation with respect to the exact 
solution. 

However, it is important to remark that, even though oq and oq are equal, the sidelobe struc¬ 
ture has a completely different physical origin. In the classical case, the sidelobes are a consequence 
of discontinuities in the surface currents, while in the quantum case the sidelobes are a result of 
quantum interference. Of course, the mathematical description of both is the same, which results in 
their similar cross sections. 

Nevertheless, these results appear to offer the prospect of exploiting the quantum sidelobe 
structure of oq to detect targets oriented off the specular direction. While current tapering will tend 
to suppress the classical sidelobes, most probably it will not have an effect on the quantum sidelobes. 
Furthermore, it is worth exploring any possible sensing mechanism that could enhance the quantum 
sidelobe structure of a target. These are subjects that we will consider in more detail in the future. 

6.5.4 TARGET’S GEOMETRY 

Let us assume now three rectangular targets of the same area, but different proportions: 

Ai = 2.5 x 4.0 X 2 (6.29) 

A 2 = 5.0 x 2.0 X 2 
At, = 10.0x1.0 A 2 

The plot of oq vs. 9 at 0 = 0 for these targets is shown in Figure 6.4. It can be seen that, even 
though they have the exact same area, the difference in their proportions alters their oq. It can be 
observed that the maximal value of oq is the same in all cases, but the sidelobe structure is different. 
This result is equivalent to the one seen for oq in Chapter 4. 

6.5.5 RANGE INDEPENDENCE 

Figure 6.5 shows that oq is indeed independent of the range to the target R, one of the most 
important features of a well-defined oq. The plot shows oq vs. 9 (4> = 0) for R/a = 1, 5, 10, 10 2 , 
10 3 , 10 4 and 10 6 . 
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Figure 6.4: Simulation of erg vs. 9 (<p — 0) for a 2D rectangular plate of the same area but different 
proportions (2.5 x 4.0 [red], 5.0 x 2.0 [green], and 10.0 x 1.0 [blue]). 

6.5.6 SMALL SIZE TARGETS 

Figure 6.6 shows erg in the case of small targets (in comparison to the radar wavelength). For a = 5 A 
we have the usual sidelobe structure. For a = A there are only two sidelobes. And for a = 0.1 k the 
quantum radar cross section is approximately constant. This behavior is also seen for the classical 
radar cross section. For small targets, the radar system cannot elucidate any details, as the target 
appears to have spherical symmetry. This final observation is not surprising, an optical-based sensor 
is only good to resolve objects larger than its wavelength. 

6.5.7 HIGH FREQUENCY PHOTONS 

The case of very high radar frequencies is shown in Figure 6.7. In the figure we show erg for a wave¬ 
length of the same size as the interatomic distance <5. The result is a plot similar to a crystallographic 
spectrum. Indeed, the quantum radar is probing the atomic structure of the target. 

Of course, this is not a desired feature of a quantum radar. This simply means that the quantum 
radar model has an upper and lower limit of valid frequencies. But then again, this is the exact same 
situation as with a classical radar. 

6.5.8 ATOMIC STRUCTURE 

Figure 6.8 explores how oq depends on the atomic structure of the target. We consider a target of the 
same dimensions, but with different number of atoms and interatomic distances. The plots depict 
erg vs. 9 (0 = 0). The cases considered are arrays of 100 x 2, 1000 x 2, and 2500 x 2 atoms with 
interatomic distances 8 = 10 —3 , 10 -4 , and 4 x 10~ 5 X, respectively. We observe a slight variation 
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Figure 6.5: oq vs. 9 {<p = 0) for R/a = I, 5, 10, 10 2 , 10 3 , 10 4 and 10 6 .The plots make evident that aQ 
is independent of the range to the target R. 



Figure 6.6: Behavior of aQ for small targets. 


in the amplitude values of oq for high incidence angles, but all the cases retain the same periodic 
structure. 
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Figure 6.7: Quantum radar cross section for high frequency photons (7. = 0.001 and 5 = 0.001). 


QRCS 



Figure 6.8: Plots of gq vs. 6 (<f> — 0) for arrays of 100 x 2, 1000 x 2, and 2500 x 2 atoms with inter¬ 
atomic distances S = 10 -3 , 10 -4 , and 4 x 10 -5 , respectively. 

6.5.9 MULTIPLE PHOTONS 

So far we have considered the structure of oq for the case of single photon signals in the quantum 
radar pulse. It is very easy to generalize the equations to describe quantum radar pulses with n y 
photons. The definition of oq is the same, but we need to obtain a new expression for {I s ( r d> 0)- 
Let us consider as an example the case of two photons, n y = 2, one with momentum p and 
the other with momentum q (and for simplicity we ignore the polarization states). We can follow 
the same steps used in Chapter 2 to obtain expressions for the probability of photodetection. In this 
case, the photodetector needs to detect two photons, one at (r, t) and another at (r\ t'). Then, the 








































































































































6.5. SIMULATION OF o Q FOR RECTANGULAR TARGETS 105 


transition amplitude of the detector signaling the detection of the two photons is given by: 

a/(r, t, r', t') = <0|£ (+) (r, t) E (+) (r', F)|l p l q ) (6.30) 

and the measured intensity is: 

(/pq(r,fr'd'» = \af(r,t,r',t')\ 2 (6.31) 

= <0|£ (+) (r,0 E^\r', t') |l p l q )| 2 

= 'FpqCr, t,r', t')\ 2 


Now, using the expressions for the quantum fields and the Fock states we get: 


vj/ pq (r, t,r',t') = 
oc 


oc 


(X 

(X 

where we have used the fact that: 


<0|£ (+) (r,0 £<+y,f')|l p lq> 

(0|ffi p e- /( “p ? -P r) + 

x (d pe -' ( ^ ? '-P' r V a q g- , 'H''-l- r ' ) )|lpl q > 

+ e -F<V-p-r) e -«'K ? '-q.r') , lplq) 

+ r iK ( - q .,) e -i(yV) (0 | 5(] . p|1 ^ ) 

+ e -iKH*) r i(«,i , -q'«') ( o | 3 q 5 q|lpli) 

g -i(a>pf-p-r) g -;(ft> 9 f , -q-r') (oj,(- q-r) e -i(6j /) r'-p-r') 

vF yp (r, t) V Y y, t') + 'Fy q (r, t) *^(r', F) 


(6.32) 


(0|d p a q | l p l q ) = (0|o q d p | l p l q ) oc 1 (6.33) 

(0|C7pC7p| 1 plq) = (0|n q C2 q | l p 1 q ) OC 0 

which is a consequence that d p (o q ) destroys a photon of momentum p (q). 

Therefore, the intensity is given by: 

</ pq (r, t, r', t')) cx |T yp (r, t) vk /q (r', t') + * y ,(r, t ) ^(r', F)| 2 (6.34) 


which is a symmetric bipartite wave function as required by Bose-Einstein statistics. Furthermore, 
notice that in general: 


</ pq (r, t,r', F)) jk |vk Kp (r,0| 2 l^ q (r',F)| 2 


(6.35) 


or equivalently: 


</pq> 7^ </p></q> 


(6.36) 


which means that the two photodetections at (r, t ) and (r', F) are not independent but are corre¬ 
lated [73]. The reason for this behavior is that each of the detection events gives some information 
about the state of the quantum field, and this information has an influence on the estimates of sub¬ 
sequent detections. As a consequence, this is a quantum mechanical effect due to the interference 
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between the photon states in the symmetric bipartite state shown in Equation (6.34). Indeed, it is 
only possible to factorize the contributing intensities for states that represent classical light. 

Let us now consider the case of three atoms (N = 3) scattering two photons ( n y = 2). In this 


case: 


(7(r rf , 0) oc - 


n yi ' n n ~ n Y2 


(6.37) 


where the wave function subindex marks the photon and the upper index identifies the scattering 
atom. That is, the two photons can interact with any of the three atoms, and all possible combinations 
contribute to the interference pattern. Notice that there is no permutation over the photon indices 
because the photons are bosons, and therefore indistinguishable from each other. 

Then, for N atoms and n y photons: 


Cuxd,t)) 


i 

7¥7 

\n Y J 


a\a2...a, 


.(«/) 

>y 1 1 Yj 

j= 1 


n*J 


(6.38) 


where a/ = 1, 2, ...N, and sum over equal “a” indices is implied. 

The simulation of multiple incoming photons has exponential complexity because of the 
binomial combinatorial summation in the equation. Indeed, the simulation of (n y + 1) photons 
requires this number of operations: 


A- \ 
n y + 1/ 


N\ 

(n y + 1)! (N-n y - 1)! 
N\ N - n y 

(n y )\ (N — n y )\ n y + 1 
/ N\ N — n y 
\n Y ) n y + 1 


(6.39) 


And, if N n y : 

« 0(N/n y ) (6.40) 

Hy + 1 

That is, the simulation of a single extra photon requires a factor of N extra operations. Then, the 
simulation of n y photons requires: 

0({N/n y ) n y) (6.41) 


operations. 

Figure 6.9 shows the plots of oq vs. 0 (0 = 0) for n y = 1,2, and 5 photons (for a rectangular 
target made of 10 x 5 atoms). It can be observed that as the number of photons increases, so does 
the peak of the maximum value of ctq at the specular orientation. That is, a target near the specular 
direction appears to look bigger when illuminated with a quantum radar. For example, we will 
observe the magnification of the quantum cross section at the specular direction when we use an 
interferometric quantum radar using NOON states. This appears to be a pure quantum-mechanical 
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Figure 6.9: Plots of oq vs. 0 (4> = 0) for n y = 1 (red), 2 (blue), and 5 (green) photons. 

effect due to the interference terms that lead to Equation (6.36). Thus, only in the case of beams of 
classical light will we expect to see the same cross sections without regard to the number of photons 
per beam. 

However, as the number of photons increases, the sidelobe structure decreases dramatically 
and the width of the specular orientation peak becomes more narrow. Such behavior is to be expected. 
Indeed, the interferometric analysis only takes into consideration specular returns and the associated 
sidelobe structure due to quantum interference. As the number of photons increases, the dynamics 
of the scattering transition to the classical domain. In such a case, we would expect that oq will 
be extremely narrow at the specular direction and without sidelobes. That is, this behavior makes 
evident the large photon number limit for oq: 

lim oq oc oc oc 8(6) (6.42) 

n y —>-00 

which corresponds to specular returns without end-region contributions. This further shows that 
the quantum side lobe structure is not related to the side lobe structure observed in the context of 
classical radar cross sections. That is, the quantum side lobe structure is a purely quantum mechanical 
feature. 

Numerical values for the peak value of the quantum radar cross section on the main (6 = 0) 
and secondary (0 = tr/4) lobes are given in Table 6.1. It almost appears that the main lobe peak 
value increases by a constant factor as a result of adding an extra photon. However, this is not the 
case. Indeed, let us define the growth ratio G r (n y ) as: 

Oq OX (n y ) 

Oq QX (n y - 1) 


G r (fty) — 


(6.43) 
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where <jQ ax (n y ) is the maximal value of gq for n y photons. Thus, we can observe a steady decrease 
in its value as we increase the number of photons. Furthermore, we can expect that the growth ratio 
will steadily decrease to zero for a large number of photons. 

On the other hand, the peak of the lateral lobes appears to decrease exponentially faster with 
the number of photons. It is also worth mentioning that the “damping” of the quantum side lobes 
due to the addition of extra photons to the radar pulse depends on the ratio N/n y . In more realistic 
examples, where N n y , the damping will not be as severe as the one shown in the figure. 


Table 6.1: Peak values of the quantum radar cross section on the main (9 = 0) 
and lateral (0 = tt/ 4) lobes. 

No. of Photons 

Main Lobe Peak 

Growth Ratio 

Lateral Lobe Peak 

1 

0.0003152 


1.12 x 10 -3 

2 

0.0006474 

2.054 

1.14 x 10 -6 

3 

0.0009444 

1.459 

8.06 x 10 -8 

4 

0.0012262 

1.298 

4.84 xl0“ 9 


At this point is worth discussing the two main quantum mechanical effects that contribute 
to the behavior of the quantum radar cross section: 

• First, the sidelobe structure is a pure quantum mechanical effect due to the interference of the 
photon wave functions and it has no classical counterpart. Indeed, classically there would be 
no sidelobes without discontinuous surface currents due to edge effects. 

• And second, the maximal value of the quantum radar cross section increases with the number 
of photons in the quantum radar beam. This effect is also due to the interference between the 
photons as expressed by Equation (6.36). 

Because of the correspondence principle, we expect that the quantum radar cross section will 
become the classical radar cross section for a large number of photons. The first effect appears to obey 
this principle, as the sidelobe structure disappears as the number of photons increases. However, the 
second effect appears to persist in the large photon number, as the maximal value of the quantum 
radar section continues to increase with the number of photons. 

At this point it is important to recall the second footnote of Section 2.5, where we mentioned 
that a large number of photons is not enough to guarantee the correspondence principle. Indeed, 
the Fock states used for the computation of the quantum radar cross section are non-classical even 
with a large number of photons. To obtain the principle of correspondence for the second effect one 
needs to explicitly use coherent states with large average photon numbers. It is only in such a case 
that the expression for the contributions to the intensity becomes an equality: 

(hi) = (h)(h) (6-44) 

which means that there are no interference effects and as a consequence there will be no change to 
the value of the radar cross section. 
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Therefore, in the small photon number regime, it appears that specular returns and quantum 
interference are the principal contributors to the value of oq . Therefore, it may be worth considering 
if it is possible to detect targets using the sidelobe structure of quantum radars operating in the small 
photon number regime. In such a case, not only an optimal value for n y will need to be found, but 
also explore strategies and protocols that amplify the quantum sidelobe structure of a target. 

6.6 SUMMARY 

In this chapter we have discussed the need to define a quantum radar cross section oq. The definition 
taken follows the same structure as the one for classical radar ac- We used interferometric arguments 
based on quantum electrodynamics to determine the intensity measured by the radar receiver. It was 
found that for single photon quantum radars, oq has a sidelobe structure which is a purely quantum 
mechanical effect (in contrast, the sidelobe structure of oq is a consequence of end-region returns). 
We also observed that oq > oq for rectangular targets near the specular direction if the quantum 
radar pulse has more than one photon. In addition, just as multiple photon pulses tend to increase 
the specular return of oq , they also narrow the peak and reduce the sidelobe structure. 

These results suggest that it may be possible to detect targets using the sidelobe structure of 
quantum radars. This is an idea that deserves to be explored in more detail. 
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CHAPTER 7 

Conclusions 


Preliminary results show that interferometric quantum radar using entangled photons can provide a 
quadratic increase of resolution over non-entangled photons even in an attenuating environment [97]. 
Quantum illumination offers higher signal-to-noise ratio which increases the detectability and iden¬ 
tification of stealth targets even in the presence of noise or a jamming signal [67]. In addition, defining 
the quantum radar cross section using quantum electrodynamics leads to the result that the “effec¬ 
tive visibility” of a target near the specular direction is increased if observed with a quantum radar 
instead of a classical radar [61, 62]. Furthermore, a quantum sidelobe structure offers the possibility 
of detecting RF stealth targets [61]. 

As a consequence, quantum radars may be able to detect, identify, and resolve RF stealth plat¬ 
forms and weapons systems. In the electronic battlefield, quantum radar may become a revolutionary 
technology just as RF stealth technology was in the last three decades of the 20th century. 

7.1 OPEN QUESTIONS 

Quantum radar is a novel technological concept and its development presents a variety of multidis¬ 
ciplinary, experimental, and theoretical research challenges. Some of these include: 

1. What is the most optimal way to entangle photons in the X-band frequency regime? 

2. What is the most optimal frequency to be used by quantum radars? 

3. What is the best physical design of a quantum radar? 

4. What is the best way to measure the operational performance of a quantum radar? 

5. How to deal with the effects of noise in the atmosphere? 

6. How to describe the capacity of a quantum radar channel? 

7. How to simulate quantum radar cross sections? 

8. What is the optimal number of coherent photons per radar pulse? 

9. How to design low observable quantum stealth platforms? 

10. How to use weak measurements to enhance detection? 

11. How to exploit other quantum resources to enhance target detection? 
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12. How to perform signal detection from quantum radar measurements? 

13. How to ruggedize a quantum radar for its deployment in an operational environment? 

14. What are the effects of geometry on the quantum radar cross section? 

15. How do these results change when using realistic targets in arbitrary motion? 

16. What type of signal analysis is needed to detect and track targets using a quantum radar? 

17. Is there a quantum sensing strategy that amplifies the quantum sidelobe structure? 

7.2 THE BOTTOM LINE 

Quantum radar is a promising technology that could have a strong impact on the civilian and military 
realms. It is important to note that, even though quantum sensing technology is not as mature as 
quantum computation, the engineering of a full-blown multi-qubit quantum computer is a much 
harder challenge than the design of deployable quantum sensors. As a consequence, quantum radar is 
a high risk, high payoff proposal. By any means, the concept of a quantum radar is worth investigating 
in further detail. 
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